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PREFACE 


This  report  presents  results  of  a  study  designed  to  Investigate 
the  dynamics  of  a  spin-stabilized  free  rocket  during  the  launch  phase 
of  motion.  The  models  developed  are  applicable  for  the  case  of  very 
high  spin  rates  similar  to  those  achieved  by  ANSSR,  an  aerodynamlcally 
neutral  sp1n-stab11l2ed  rocket  developed  by  Emerson  Electric  Company 
for  the  U.  S.  Anqy  [lj*.  This  105  mm  rocket  Is  pre-spun  In  Its  launcher 
by  four  auxiliary  motors  thrusting  tangentially  at  the  rear  of  the  round. 
The  boost  phase  of  the  flight  does  not  begin  until  these  motors  have 
been  exhausted.  Thus,  the  rocket  exits  the  launcher  with  spin  rates 
(10,000  rpm)  considerably  greater  than  those  considered  by  previous 
Investigators  [2].  A  double-diameter  launcher  tube  Is  used  to  minimize 
tip-off  effects. 

The  design  objective  for  the  ANSSR  system  his  been  to  utilize  the 
stabilizing  effect  of  the  extremely  high  spin  rate  in  lieu  of  fins  to 
decrease  the  overall  wind  sensitivity  of  the  weapon.  A  factor  affecting 
the  accuracy  of  any  free  rocket  Is  mallaunch  (that  component  of  the  rocket 
angular  rate  vector  not  aligned  with  the  launcher  bore  axis  at  the  instant 
of  release).  Mallaunch  values  observed  with  optical  lever  and  camera 
techniques  at  Redstone  Arsenal  ware  higher  than  expected  [3],  This 
investigation  has  addressed  the  problem  of  mallaunch  due  to 


*  Numbers  in  square  brackets  refer  to  references  at  the  end  of  the  preface. 
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•  tip-off  effects 


•  In-tube  rocket  bending 

•  launcher  motion  caused  by  rocket 

excitation 

This  report  Is  presented  In  three  distinct  sections.  Each 
section  Is  complete  with  Its  own  list  of  symbols,  fibres,  tables 
and  references.  Results  obtained  are  discussed  at  the  end  of  each 
section  and  will  not  be  belabored  here.  However,  we  will  mention 
that  the  results  show 


•  that  tip-off  should  not  be  a  major 
ca'ise  of  mal launch 


•  that  body  bending  can  be  a  major 
cause  of  mal launch 


•  that  the  launcher  support  system 
can  be  designed  so  as  to  minimize 
mal launch 

The  diligent  efforts  of  the  project  leader's  co-workers  Messrs. 
Lindholm  and  Tomb  is  hereby  acknowledged.  The  calculations  were 
made  with  the  help  of  Mr.  John  B.  Hill.  Also,  the  skill  and  extreme 
good  humor  of  Mrs.  Dianne  Fretwell  along  with  Mr.  Robert  Sullivan 
and  Mr.  Robert  Young  were  indispensable  In  preparing  the  manuscript. 

The  research  team,  during  the  course  of  this  Investigation,  has 
enjoyed  the  cooperation  of  the  Aeroballlstlc  Directorate  of  Redstone 
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Arsenti  and  particularly  that  of  Messrs.  Conard  and  Vest  of  the  Systems 
Evaluation  Branch. 


Joe  W.  Reece 
Project  Leader 
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SECTION  I 


MALLAUNCH  DUE  TO  TIP-OFF 

ABSTRACT 

A  lineanized  model.  is  developed  {a*  die  dynamic s  o{  a 
spinning,  unbalanced,  ligid  locket  doling  the  tip-o {{  phase 
o{  launch.  Mallaunch  iaX.es  {,01  the  locket  aie  obtained  {on 
the  shoot  time  int envoi  {oom  tip-o {{  to  end  o{  guidance.  The 
oesults  show  that  tip-o {{  should  not  be  a  majol  cause  o{ 
mallaunch  {of i  veoy  hi^h  spin  nates  and  shont  tip-o{{  times. 
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LIST  OF  SYMBOLS 


X.Y.Z 

U.k 
x.y,  2 

*1 ,®2’*3 


^x^yy^zz 


ground  reference 
unit  vectors  of  X.Y.Z 
body  reference 

unit  vectors  of  x,y,z 

Euler  Angles 

angular  rotation  vector  of  x,y.z  with  respect  to  X.Y.Z 

taH;ya:zr°tati0n  V6Ct0r  °f  the  rocket  ^pect 

angular  rotation  vector  of  the  rocket  with  result 
to  x,y,z  (spin  rate  of  rocket)  respect 

xv:,rtr of  *• "*•* 

moment  re0401*  "1th  respect  to  the  origin 

raJIatoTy??  0f  inert1a  of  the  roeket  w1th 

front  bearing  reaction  vector 
thrust  vector 
inertia  load  vector 
rocket  mass 

distance  of  rocket  mass  center  from  y-axis 
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position  vector  of  the  c.m.  of  the  rocket  with 
respect  to  x,y,z 

components  of  ff0  expressed  In  x,y.z 
32.2  ft/sec2 

transfer  functions  of©  and  V 
eigenvalues  of  the  system 

magnitude  of  the  pitch  mal launch  rate 

magnitude  of  the  yaw  mal launch  rate 

approximate  rocket  velocity  at  the  end  of  guidance 

average  rocket  acceleration  from  boost  motor  Ignition 
until  the  end  of  guidance 

guidance  length 

bearing  mismatch  -  longitudinal 
total  tip-off  time 


A.  INTRODUCTION 


A  possible  source  of  mal launch  for  the  ANSSR  Rocket  Is  due  to  what 
has  come  to  be  called  the  "tip-off"  effect.  This  phenomenon  occurs  In 
the  time  period  after  the  rocket's  rear  bearing  has  lost  Its  launcher 
tube  support  while  the  front  bearing  retains  Its  support.  Due  to  the 
extremely  small  time  of  contact  It  Is  possible  to  obtain  a  closed  form 
solution  for  a  linearized  dynamical  model  of  the  rocket  motion  during 
"tip-off".  The  process  of  analysis  allows  the  capability  of  treating 
both  static  and  dynamic  unbalances.  The  mallaunch  rates,  caused  solely 
by  tip-off,  which  are  predicted  by  this  model  are  the  pitch  and  yaw 
rotation  rates  for  the  rocket  Immediately  after  the  front  bearing  has 
lost  Its  launch  tube  support,  i.e. ,  as  the  rocket  becomes  free.  These 
rates  are  depicted  In  Figure  1. 

The  spinning,  unbalanced,  rigid  (I.e.,  no  bending)  rocket  during 
the  tip-off  phase  Is  open  to  the  same  sort  of  dynamic  analysis  usually 
applied  to  the  classical  dropped  gyro  [1],  A  significant  difference 
Is  the  inclusion  of  the  pertubation  of  the  rocket  motion  due  to  static 
and  dynamic  mass  unbalances.  As  a  result,  this  analysis  yields  a  set 
of  linearized  equations  of  motion  rather  than  the  more  general  solution 
utilizing  elliptic  functions  [1]. 

Shown  in  Figure  2  is  the  configuration  of  the  rocket  in  the  tip-off 
phase.  In  the  analysis,  the  rocket  is  assured  to  be  simply  supported 
by  its  front  bearing.  These  simple  supports  are  aligned  with  the 
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gure 
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geometric  center  of  the  front  bearing.  At  "time  zero",  (l.e. ,  when  the 
rear  bearing  has  lost  Its  launcher  tube  support)  the  position  of  the 
rocket  Is  as  shown  In  Figure  2. 

The  physical  characteristics  used  In  this  analysis  are  presented 
In  Table  I. 


Rocket  Mass  (Sluq) 

1 

Ensnssi 

1.54 

2.2 

0.03 

2.058 

TABLE  I.  Rocket  Characteristics  (ANSSR  II) 


B.  ANALYSIS 


The  position  of  the  rocket  In  space  can,  at  any  time,  be  described 
by  the  two  Independent  Euler  angles  shown  In  Figure  3. 

The  coordinate  systems  and  Euler  angles  used  In  the  analysis  are 
defined  as  follows: 

•  Ground  Reference:  Ground  reference  Is  (X,Y,Z)  with  respective 
unit  vectors  (i.j.k).  The  origin  Is  fixed  at  the  geometric  center  of 
the  front  bearing.  The  Y  axis  Is  taken  positive  going  out  of  the 
launcher  tube's  longitudinal  axis.  The  Z  axis  Is  taken  positive  up  and 
parallel  to  the  local  vertical.  The  X  axis  completes  a  right  hand 
system.  All  equations  of  motion  will  be  written  with  respect  to  this 
coordinate  system.  This  can  be  done  as  the  acceleration  of  point  "o" 

Is  essentially  parallel  to  the  acceleration  of  the  rocket  C.M.  for 

the  short  time  period  Involved.  This  frame  will  be  designated  as 
"frame  g". 

•  Body  Reference:  Body  reference  Is  (x,y,z)  with  respective  unit 
vectors  (e^,e2,e3).  The  origin  (point  "o")  of  this  frame  Is  fixed 
coincident  with  that  of  frame  g.  At  time  zero  this  frame  Is  aligned 
with  frame  g  (l.e.  ,♦*  9  *  0.0).  At  any  time,  the  position  of  the 
body  reference  with  respect  to  frame  g  is  defined  by  the  Euler  angles 

♦  and  9  as  In  Figure  3.  This  frame  will  be  designated  as  "frane  b". 

It  Is  Important  to  note  that  x  and  z  do  not  spin  with  the  rocket, 
l.e. ,  the  rocket  spins  about  the  y  axis  of  frame  b. 
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•  Systam  Euler  Angles:  The  two  Independent  Euler  angles  that 
describe  the  position  of  frame  b  with  respect  to  frame  g  are  #  and  #. 
♦  Is  the  first  angle  of  rotation  about  Z  and  0  Is  the  second  angle  of 
rotation  about  x.  Both  angles  and  their  rates  are  expressed  In  rad 
and  rad/sec  respectively. 

The  relation  between  the  coordinates  of  frame  b  and  those  of 
frame  g  Is: 


with  matrix  [A]  defined  as: 

'10  0*1  ft*  S*  O' 

[A]  -  0  CO  SO  •  -SO  C#  0  (2) 

0  -SO  CO  Lo  0  1. 


where  CO  Is  taken  to  be  COS(O),  etc. 

Now: 

g— R  -  the  angular  rate  of  frame  b  with  respect  to  frame 
g  (m  rad/sec). 

g^R  •  the  angular  rate  of  the  rocket  with  respect  to  frame 
g  (m  rad/sec). 

O  -  the  angular  rate  of  the  rocket  with  respect  to  frame 
b  (f  rad/sec). 

The  above  angular  rates  are  given  by  the  following  equations: 


g— b  ■  0  ej  +0S0  eg  +OC0  e.j 


(3) 
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9S*  "  •  *1  *  <♦  S#  +  n  )  «2  +^C#  *3  (*) 

0  •  «  e2  (5) 

It  follows  that  the  angular  momentum  vector  of  the  rocket  expressed 
In  frae  b  Is: 

■  '«*  *1  ‘  'yy  <*S*  ♦«)  h  *  '«♦<:#  ?3  (6) 

The  sun  of  all  external  moments  about  point  *0"  must  be  equal  to 
the  first  time  derivative  of  HQ  with  respect  to  frame  g.  So: 


(7) 

(8) 


If  It  Is  assumed  that  for  the  time  period  of  Interest,#,  ♦  and 
their  respective  time  derivatives  are  very  small  (on  the  order  of 
ml  111  radians),  all  second  order  terms,  such  as  l^+StCI  ,  that  result 
from  the  evaluation  of  equation  (8)  can  be  Ignored  as  |«|  Is  large. 

If  S#*  9  and  C#  ■  1,  equation  (8)  reduces  to 


(9) 
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These  equations  are  similar  In  font  to  those  presented  In  [23  for 
the  m>  a  lysis  of  the  two-axis  gyroscope. 

The  only  forces  acting  on  the  rocket  during  the  tip-off  configur¬ 
ation,  Figure  4,  are  rocket  weight.  Inertia  loads  due  to  static  and 
dynamic  mass  unbalances,  thrust,  and  front  bearing  reactions.  Of  these, 
only  the  rocket  weight  and  Inertia  loads  exert  any  moment  about  point 
"0".  This  Is  true  because  the  rocket  Is  assumed  to  be  rigid  (l.e. ,  no 
bending)  and  point  "0"  Is  located  at  the  geometric  center  of  the  front 
bearing.  It  Is  known,  that  If  the  rocket  Is  bent  during  tip-off,  that 
a  component  of  the  thrust  vector  will  create  a  moment  about  point  "0”. 
The  question  of  the  degree  of  bending  that  the  rocket  experiences  Is 
treated  In  Section  II  of  this  report. 

IT  -  front  bearing  reaction.  This  force  always  acts  In  the  z  -  x 
plane  and  Is  assumed  to  act,  at  all  time,  on  a  line  through 
point  -0"  ($  lb). 

T  ■  rocket  thrust  vector.  This  force  always  acts  along  y.  lb), 

-mgk  ■  rocket  weight.  This  force  always  acts  at  the  rocket  c.m. 

and  Is  always  parallel  to  Z.  (J*  lb).  Note,  g  -  32.2  ft/sec2. 

Fj-  Inertia  load  due  to  static  and  dynamic  mass  unbalances.  This 
force  Is  assumed  to  act  through  the  rocket  c.m.  and  always  In 
a  plane  parallel  to  the  z-x  plane.  The  rate  of  rotation  of 

^1  is  assumed  to  be  that  of  the  rocket  (l.e.,  987  rad/sec) 
and  Its  location  with  respect  to  an  arbitrary  reference  axis 
In  Its  plane  of  rotation  will  be  given  by  the  quantity  (It 
where  fi*  987  rad/sec  and  t  *  tlm  (2  sec).  The  magnitude  of 
Is  given  by: 
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tj  -  mgr  -  «£  fl2  r  cos(Qt) 

04) 

-  I  #  •  «E  Q2  r  sln(Ot) 

Since  the  rocket  Is  symmetric  (l.e. ,  ■  I22)  the  following 

terms  are  defined: 


M 


1 


(15) 


H 


2 


HEQ2r 


Note  that  the  term  m  Is  the  natural  frequency  of  the  system 
(if  rad/sec).  If  equations  (15)  are  substituted  Into  equations  (14) 
then  the  equations  of  motion  take  on  the  following  form: 


•  -«+  -  m1  -  m2  cos(nt) 
♦+**#  -  m2  sln(Qt) 


(16) 


The  equations  of  motion  for  the  system  are  a  pair  of  second  order, 
coupled,  linear,  non-homogeneous ,  ordinary  differential  equations  In# 
and  +  . 

Equations  (16)  can  be  modified  for  the  case  of  Fj  being  Initially 
aligned  with  the  -Z  axis  at  time  zero  and  have  the  form: 


•  -m+  -  M1  +  M2  cos  (a  t) 
^  +•*#  *  -«2  s1n( fl  t) 


(17) 
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The  Laplace  transform  method  [2]  was  used  to  solve  equations  (16). 
The  transformed  equations  have  the  form: 


M1  M2S 
SUQ* 

h2o 

■  S2+Q2 


(18) 


where  the  S's  are  the  eigenvalues  of  (16). 

The  response  functions  for  8  and  V  can  be  obtained  from  (18) 
and  are  of  the  form: 


6 


’1 


S(S2  +  m2) 


sm2  m2o« 

(S2  +  a»2)(S2  +  n2)  S(S2  +  «2)(S2  +  «2) 


M2(fl+**)  H^e» 

*  (S2  +  «2)(S2  +  fl2)  S2(S2  +*»2) 


(19) 


By  use  of  the  partial  fraction  method  and  tables  In  [2]  the  Inverse 
transforms  of  0  and  ¥  are  found  and  the  time  responses  for  a  0  and  V 
during  tip-off  are: 


,(t)3C0S(nt)r-J2 — ].cos(*t)r!!i  + J!? — +  i 

Ln(n-*)J  L**2  «(n-«)J  *2  “Q 

♦  (t)  =  s1n(  m  t)T—  + - - - 1  -  sin  (n  t)  [ - - - 1  -  1 

l*2  «*(n -e»)J  Ln(a-e»)J  “ 

With  the  first  time  derivatives  of  equations  (19)  and  the  projection 


of  these  quantities  onto  frame  g  the  magnitudes  of  the  pitch  and  yaw 
mal launch  rates  are  found  to  be: 
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jeos(?) 

f 


(21) 


Figure  9  wes  obtained  by  assuring  that  both  the  thrust  and  the 
■ass  are  constant  from  boost  rotor  Ignition  until  the  end  of  guidance. 
An  average  thrust  for  this  tine  was  conputed  as  4808.367  1b  [4].  With 
a  rocket  nass  of  1.S4  slug  an  average  acceleration  of  3122.316  ft/sec2 
can  be  calculated.  As  the  Initial  velocity  of  the  rocket  Is  zero  the 
velocity  at  the  end  of  guidance  can  be  calculated  from: 


V  -  yflas  (22) 

and  has  a  value  of  1246.96  in/sec. 

a  ■  average  acceleration  up  to  the  end  of  guidance  (£  ft/sec2) 
s  »  guidance  length  (20.75  in.)  [4] 

If  the  velocity  of  the  rocket  Is  assumed  to  be  constant  during  tip-off 
then  the  following  relation  exists  between  velocity,  bearing  mismatch, 
and  total  t  me  of  tip-off: 

t  -  8/V  (23) 


where 

8»  bearing  mismatch  (£  In.) 
t  *  total  tip-off  time  (»  sec. ) 

V  *  average  velocity  of  the  rocket  during  tip-off  (1246.96  In/sec) 

Thus,  for  any  specified  amount  of  bearing  mismatch  the  total  tip-off 
time  can  be  calculated. 
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A  typical  dynamic  unbalance  at  the  front  bearing  (before  balancing) 
has  been  measured  as  0.45  oz-ln  and  reported  In  [3].  This  Is  equivalent 
to  an  effect  caused  by  a  D'Alembert's  force  at  the  center  of  mass  and 
expressed  by  equation  (10).  For  a  rocket  mass  of  1.54  slugs  and  (1  of 
987  rad/sec.  the  equivalent  E  trill  be  5.015  x  10"®ft. 

In  order  to  obtain  the  mallaunch  rates,  the  amount  of  bearing  mis¬ 
match  must  be  specified.  If  the  bearing  mismatch  Is  as  much  as  a  tenth 
of  an  Inch  then  the  total  tip-off  time  Is  8.0  x  10"^  sec  (Figure  9). 

With  the  aforementioned  conditions  and  for  the  case  where  Tj  Is 
Initially  aligned  with  -Z  the  pitch  and  yaw  mallaunch  rates  of 
0.048  x  10"2  rad/sec  and  -0.0136  x  IQ"3  rad/sec  respectively  (Figures  7 
and  8).  These  two  values  yield  by  a  vector  sum  a  total  mallaunch  rate 
of  0.48  x  10"3  rad/sec  due  to  the  .10  inch  mismatch. 

The  shaded  area  on  Figures  5,  6,  7,  and  8  depict  the  range  of 
mallaunch  for  0.0  to  0.02  oz-ln  of  mass  unbalance  as  reported  In  [3]. 

For  a  bearing  mismatch  of  0.23  in  and  total  tip-off  time  of  20  x  10‘® 
sec,  a  mass  unbalance  of  0.02  oz-ln  yields  pitch  and  yaw  mallaunch  rates 
of  0.0598  x  10  2  rad/sec  and  0.0123  x  1G~3  rad/sec  respectively  If  Tj  Is 
assumed  to  be  initially  aligned  with  the  -Z  axis.  These  two  values 
yield  by  a  vector  sum  a  total  mallaunch  rate  of  0.598  x  10  rad/sec 
due  ♦’c  the  0.23  inch  mismatch. 


Yaw  Hal  launch  Rate  (RAO/SCC)  x  10 
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SECTION  II 


HALLAUNCH  DUE  TO  IN- TUBE  ROCKET  BENDING 

ABSTRACT 

Rochet  {legibility  is  modelled  {on  the  launch  phase.  The 
model  account*  {on  beaning  cleanances,  a  vaniabte  sti{{ness 
mo  ton  lection  and  allow  {on  the  3-0  motion  o(  a  nigid  u>an- 
head.  The  phenomenon  o{  uhinling  is  investigated  and  the 
bent  con{igunation  o{  the  nocket  is  calculated  {on  synchnonous 
wheeling.  The  nesulting  mallaunch  nates  agnee  utith  neponted 
nange  data. 
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LIST  Of  SYMBOLS 


«  acceleration  of  payload  center  of  Mass  with  respect  to  X,Y,Z 

longitudinal  (along  launch  tube)  acceleration  of  rocket  due 
to  engine  thrust 

e^.e^e^  unit  vectors  In  x,y,z 

ff  aonent  of  nomentua  vector 

1,j,k  unit  vectors  In  X.Y.Z 

I^j  Inertia  tensor  of  rocket 

k  dynamic  stiffness 

K  stiffness  (as  to  bending) 

L  length  between  bearing  centers 
m  mass 

ff  bendl ng  moment 


r  position  vector  of  center  mass  of  payload  section  with  respect 
to  X.Y.Z 

x,y,z  body  reference 

X.Y.Z  ground  reference 

a  bending  angle 

rocket  section  deviation  angle  for  1th  section 
*^,*2  bearing  clearances 

$  angle  between  longitudinal  axes  of  launch  tube  and  unbent  rocket 

¥,e  Euler  Angles 

Mr  Laplace  transform  of  * 

9  Laplace  transform  of  e 
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«  body  fro at  spin  rat*  with  rtspoct  to  X»Y,Z 

Hallaunch  or  ccoponont  of  rockot  spin  vtctor  not  along  launch 
tuba  axis 

angular  rat*  vtctor  of  rocket  with  rtsptct  to  X,Y,Z 
o  angular  rat*  vtctor  of  rocktt  with  rtsptct  to  x,y, z 


A.  INTRODUCTION 


The  objective  of  this  section  of  the  report  Is  to  develop  an  analyti¬ 
cal  model  that  will  predict  the  whirling  and  mal launch  rates  of  ANSSR  II. 
The  motor  section  of  the  rocket  will  be  modelled  as  a  variable  El  bean 
that  Is  simply  supported  between  Its  front  and  rear  bearings  while  the 
payload  section  will  be  considered  rigid.  The  basis  of  these  asswptlons 
Is  Figure  25  of  [1]. 

The  analysis  will  be  divided  Into  four  parts.  The  bending  and 
deflection  equations  of  the  motor  section  subjected  to  the  loading  condi¬ 
tion  shown  In  Figure  1  will  be  developed  first.  A  general  motion  analysis 
of  the  payload  section,  coupled  with  the  bending  equations,  will  then  be 
carried  out  and  will  yield  a  whirling  rate  for  ANSSR  II  that  Is  within  18 
to  20  percent  of  the  observed  whirling  rate.  A  prescribed  motion  analysis 
of  the  payload  section,  taking  Into  account  the  rocket's  bearing  clearances, 
is  also  coupled  with  the  bending  equations  yielding  the  degree  of  "In-tube" 
bending  for  ANSSR  II.  The  degree  of  bending  predicted  by  this  model  will 
be  in  agreement  with  bending  data  that  is  presented  In  Figure  22  and  Figure 
23  of  [1].  Once  the  degree  of  "in-tube"  bendlno  is  known  then  the  mass 
center  of  the  bent  rocket  and  its  PLAI  with  respect  to  its  mass  center  are 
both  found  via  a  lumped  parameter  model.  Thus,  knowing  the  bent  rocket's 
PLAI,  the  mallaunch  rate  of  ANSSR  II  is  easily  found. 

The  majority  of  the  above  analysis  is  carried  out  In  one  of  the  trans¬ 
verse  planes  of  the  rocket  as  complete  synmetry  of  rocket  motion  In  its 
perpendicular  transverse  planes  is  assumed. 


28 


B.  BENDING  MO  DEFLECTION  EQUATIONS  OF  THE  MOTOR  SECTION 


The  aotor  section  of  ANSSR  II  Is  considered  to  be  •  verleble  El  bean 
as  shown  In  Figure  1.  The  coordlnete  system  used  In  this  enelysls  Is 
(7,Y,7).  The  origin  Is  located  at  the  center  of  the  rear  bearing.  The  7 
axis  Is  taken  positive  going  froa  the  rear  of  the  rocket  towards  the  nose 
and  Is  parallel  to  the  unbent  rocket's  longitudinal  axis.  The  Y  and  7 
axes  are  noraal  to  7  and  fora  a  light  hand  systaa.  The  sign  convention 
that  will  be  used  Is  that  a  negative  moment,  applied  as  In  Figure  1,  will 
cause  a  positive  displacement. 

The  bending  and  deflection  equations  for  the  various  sections  of  the 
■otor  section  are  represented  as; 


M_ 

x 

7TI 


1 


-2 

y 


ci 

^7 


(1) 


and 

Mx  _3  ci  _ 

*i  *  •  6  eii  u y  +  117  y  *  117  (2) 

where. 


1*1,2,  and  3  and  represents  the  beam  sections: 


0  iy  h*  H  iy  1 4’  4  -y  - L 

respectl vely, 

M_  *  A  noment  applied  In  the  negative  X  direction  as  shown  In 

x  , 

Figure  1  (2  ft  lb). 
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EI^  ■  the  El  of  the  1th  section  of  the  no  tor  section,  see  Figure  1 
(2  lb  ft*), 

L  ■  the  length  of  the  eotor  section  (J  ft) 

C1  «n<1  Ki  •  constants  of  Integration  of  the  1th  eotor  section  (J  lb  ft2  and 
lb  ft  respectively). 

The  constants  Cj  and  Kj  can  be  evaluated  due  to  the  following  boundary 

conditions: 

At  y  ■  0,  *  •  0; 

At  y  -  L,,  o  •  o_  ,  and  «,  » 

xT  x2  •  z 

At  y  ■  L-,  a  ■  a  ,  and  c.  ■  a,; 
z  x2  x3  2  J 

At  y  -  L,  *3  ■  0 


Upon  evaluation  of  the  boundary  conditions  the  Cj  and  K.  are' found  to  be: 
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The  degree  of  bending  In  both  transverse  planes  of  the  motor  section 
at  the  front  bearing  can  be  found  by  combining  equations  (1)  and  (5). 

The  equations  are  of  the  following  form; 


a 

X3 


K  X 


-  1  M 

a  *  -  7 

13  if  Z 


(9) 

00) 


where 


i. 

if 


OD 


The  quantity  K  can  be  thought  of  as  the  bending  stiffness  of  the 
motor  section  at  the  front  bearing  to  a  moment  applied  as  shown  in  Figure  1. 
The  numerical  values  of  lj,  Lg.  L  and  EI^  used  In  this  report  are 


presented  In  Table  I. 


C.  GENERAL  MOTION  ANALYSIS  OF  THE  PAYLOAD  SECTION 


As  stated,  the  payload  section  of  the  rocket  will  be  assumed  to  be 
rigid.  The  position  of  the  payload  section  In  the  space  can  be  described 
by  two  Independent  Euler  angles  as  shown  In  Figure  2. 

The  coordinate  systems  and  Euler  angles  used  In  this  portion  of  the 
analysis  are  defined  as  follows: 

e  Ground  Reference:  Ground  reference  Is  (X,Y,Z)  with  respective 
unit  vectors  (1,J,k).  The  origin  Is  fixed  at  the  center  of  the  front 
bearing.  The  Y  axis  Is  taken  positive  going  out  of  the  launch  tube 
longitudinal  axis.  The  Z  axis  Is  taken  positive  up  and  parallel  to  the 
local  vertical.  The  X  axis  completes  the  right  hand  system. 

e  Body  Reference:  Body  reference  Is  (x,y,z)  with  respective  unit 
vectors  (e^,  ?£,  e3).  The  origin  (CMp)  Is  located  at  the  mass  center  of 
the  payload  section.  If  the  angles  +  and  0  were  both  zero  then  the  body 
reference  axes  would  be  parallel  to  the  corresponding  ground  reference  axes. 

It  Is  important  to  note  that  x  and  z  do  not  spin  with  the  rocket, 
l.e.,  the  rocket  spins  about  the  y  axis  of  the  body  frame. 

•  System  Euler  Angles:  The  two  Independent  Euler  angles  that  des- 
scrlbe  the  position  of  the  body  reference  (or  frame)  with  respect  to  the 
ground  reference  are  ♦  and  0.-#  Is  the  first  rotation  about  z  and  -0 
Is  the  second  angle  of  rotation  about  x'.  Note  that  x'  Is  parallel  to  x. 
Both  angles  and  their  time  derivatives  are  expressed  In  rad  and  rad  per 
second  respectively. 
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The  relation  between  the  coordinates  of  the  body  frmm  and  those  of 
the  ground  reference  Is: 

[i]  •"’•[;] 

with  the  matrix  [8]  defined  as. 


or,  the  relation  may  also  have  the  following  form: 


03) 


04) 


0  0  ’ 

C»  S#  (15) 

-S#  C# . 


where  C#  Is  taken  to  be  COS(i),  etc. 

Now: 

u  *  the  angular  rate  vector  of  the  body  frame  with  respect  to  the 
ground  reference  (jj  rad/sec). 
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^  ■  tha  angular  rata  victor  of  the  rocket  with  respact  to  the  ground 
reference  (2  rad/sec). 

Q  »  the  angular  rate  vector  of  the  rocket  with  respect  to  the  body 
frame,  l.e. ,  the  rocket  spin  rate  (4  rad/sec). 

The  above  angular  rate  vectors  are  given  by  the  following  equations: 


id  *  -  de^ 

+  d>S#  e2  -  d*  C#e3 

06) 

N 

1 

J?f 

♦  (dSd  +  ft)e2  -  ♦Cdi’j 

07) 

ft  *  ne2 

(18) 

It  follows  that  the  angular  momentun  vector  of  the  payload  section 

with  respect  to  Its  mass  center,  expressed  In  the  body  frame  Is  [2]: 

ffanp  ■  ■  *  'yy  <*»  *  “)'2  -  !ZI  * «  <19> 

The  sum  of  all  external  moments  about  the  mass  center  of  the  payload 

section  must  be  equal  to  the  first  time  derivative  of  with  respect  to 

the  ground  reference.  Therefore, 

■  *[<%>]  Body  *"*%  <20> 

frame 

A  constant  spin  rate  of  987  rad/sec  will  be  used  in  this  analysis. 

This  value  was  obtained  by  averaaing  the  spin  rates  (excluding  flights 

1,  2,  and  3)  presented  In  Table  6  of  [3].  As  987  rad/sec  is  a  very  large 

•  2 

number,  all  second  order  terms,  such  as  I2Z  d1  SB  C#  ,  that  result  from 
the  evaluation  of  equation  (20)  can  be  ignored.  Because  the  degree  of 

/ 
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"in- tube"  bending  has  been  observed  to  be  very  small,  [1  or  2].  the  sub¬ 
stitution  of  S#»  9,  S#>  C#  »  1,  and  C+*  1  Into  equation  (20)  yields: 


*o*>  "  (Iyy  a+’  *yy  Q  •  >*3  {21) 

where 

*xx  *  *zz  "  ^  transverse  mass  moments  of  Inertia  of  the  payload 
section  with  respect  to  Its  mass  center  (•  slug  ft2). 

I  ■  the  longitudinal  mass  moment  of  Inertia  of  the  payload 
section  with  respect  to  Its  mass  center  (£  slug  ft2). 
Equation  (21)  can  be  expressed  In  the  ground  reference  with  the  use  of 
equations  (14)  and  (15).  Therefore, 

%,  ■  11  yy  “  *  '  ’*,*>'  ‘  <'»♦*  'yy  «  **  <22> 

The  only  external  forces  that  create  a  moment  about  the  mass  center  of 
the  payload  section  are  the  reaction  forces  acting  at  point  "0"  as  shown 
in  Figure  3.  These  forces  can  be  evaluated  with  the  use  of  Newton's 
Second  Law,  l.e. , 


"p1 


(23) 


where 

M  *  the  mass  of  the  payload  section  (£  slug). 

1  •  the  acceleration  of  the  mass  center  of  the  payload  section 
with  respect  to  the  ground  reference  (jj  ft/sec2). 

Now,  a  can  be  found  due  to  the  fact  that. 


a 


(24) 
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where 

r  ■  the  position  vector  of  the  mass  center  of  the  payload  section 
with  respect  to  the  ground  reference  (2  ft)  and  Is  of  the  form, 

r  -  t+1  ♦  i  j  -  t«k  (25) 

where 

t  •  the  linear  distance  between  the  point  “OTp*  and  the  point  "0", 
see  Figure  2  ft).  If  equation  (25)  Is  combined  with  equation  (24) 
then  the  following  expression  Is  found, 

•  (26) 

As  the  rocket  Is  also  accelerating  In  the  positive  Y  direction  due  to  the 
thrust  provided  by  the  main  engine,  equation  (26)  can  be  modified  to  be 
of  the  form, 

a  *  t?1  +  aT  j  -  iS k  (27) 

where 

aj  *  the  longitudinal  acceleration  of  the  rocket  due  to  the  thrust  of 
the  main  engine  (*  ft/sec2). 

The  mass  of  the  rocket  and  the  thrust  are  assumed  constant  from  boost 
motor  ignition  until  the  end  of  guidance.  An  average  thrust  for  this  time 
was  computed  as  4808  lb  [4].  With  a  rocket  mass  of  1.54  slug  [5]  an 
average  acceleration  (a^)  of  3122  ft/sec2  can  be  calculated. 

The  only  external  moments  that  act  on  the  payload  section  are  the 
reaction  moments  acting  at  point  "0M,  see  Figure  3.  These  moments  can 
be  evaluated  with  the  use  of  the  following  equation, 

R  '  an  on 
P  P 


(28) 
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where 

Hg  •  the  reaction  moments  acting  on  the  payload  section  at  point 
-0-  (d  ft  1b).  . 

An  expression  for  ffR  can  be  found  by  substituting  equations  (22) ,  (23), 
(25),  and  (27)  Into  equation  (28).  This  substitution  yields  a  vector 
equation  for  Rr  which  can  be  represented  by  the  following  scalar  equations: 


HPX  *  12  V  "p  *T  *♦»  «yy 


"Ri-*(Ixx,,2V**tMP*T-»0Iyy  l31) 

The  above  equations  are  the  components  of  the  reaction  moment,  acting 
on  the  payload  section  at  point  M0N,  that  enables  the  payload  section  to 
move  as  depicted  In  Figure  2.  An  equal,  but  opposite,  moment  must  act 
on  the  motor  section  and  is  of  the  form. 


MX  =  *(Ixx  +  *2  «p)  Mp  aT  -#o  Iyy  (32) 

«Y  3  0  (33) 

MZ  3  *(Ixx  +  Hp)  -  M?  aT  +  in  Iyy  (34) 

If  the  unbent  rocket  Is  placed  in  the  tube  such  that  the  (XYy.T)  axes 
are  parallel  to  their  respective  (X,Y,Z)  axes,  and  equations  (9)  and  (10) 
are  re-written  as: 


M  *  K  a 
I  x3 


(9a) 
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y 

then  the  following  expressions  result: 

K  "j,  ■  *(•„  *  *2  v  -  •« ",  ‘ 

K  “jj  ■  '^’xx  *  *2  V  -  ♦*  %  *1 


no*) 


K'xx  ♦  *2  V  •  **  \  *T -♦«  ‘yy 

(35) 

ff'xx  *  ‘2  V  •  ♦*  MP  *T  4  •»  'yy 

(36) 

As  #  and  #  are  a  and  a  ,  to  the  first  order,  respectively, 

X3  13 

equations  (35)  and  (36)  become 

*(Ixx  +  12  mp)  '  i(l  mp  aT  +  “  0  (35a) 


mn  *  *2  "ip)  -  *  (i  »p  ^  +  K)  +  9  a  I  -  0  (36a) 

The  term  t  aip  aT  <<  K  and  can  be  Ignored  In  this  analysis.  If 

*  ■  - - — -  (3/) 

‘xx  4  "2  "p 


I  ♦  t2  m 
xx  p 


then  equations  (35a)  and  (36a)  become, 

•  •  A#  -  B i  *  0 
and 

#  -  A*  +  8  0  *  0 


(35b) 


(36b) 
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The  above  procedure.  I.e. ,  the  combination  of  general  dyn aeries  equa¬ 
tions  and  bending  equations  Is  similar  to  a  procedure  that  Is  presented 
In  [6]. 

The  Laplace  Transform  Method  [7]  can  be  used  to  put  equations  (35b) 
and  (36b)  Into  the  following  form: 

'(S2  -  A)  -  BS  *1  fel 

L  -  [C]  (39) 

BS  (S2  -  A)  J  |vj 


where 

S  *  the  eigenvalues  of  (35)  and  (36)  (•  rad/sec). 

0  and  ¥  *  the  transfer  functions  of  f  and  #  respectively. 

[C]  *  A  matrix  that  Is  a  function  of  S  and  the  Initial  values  of  #, 

•  • 

*,  t,  and  a  [7] -  Its  particular  form  Is  of  no  Importance  to 
this  analysis. 

The  eigenvalues,  or  natural  frequencies  of  whirl,  of  equations  (35)  and 
(36)  can  be  found  by  setting  the  determinant  of  the  matrix 


(S2  -  A)  -  BS 
BS  (S2  -  A)  _ 


(40) 


equal  to  zero.  This  yields  the  characteristic  equation  of  equations  (35) 
and  (36)  whose  roots  are  the  eigenvalues.  The  characteristic  equation  Is 


(S2  -  A)2  +  B2  S2  *  0  (41) 

The  physical  properties  of  the  motor  section  and  the  payload 
section  (Tables  I  and  II)and  the  spin  rate  of  the  rocket  (987  rod/sec) 
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are  all  that  Is  needed  to  evaluate  equation  (41).  The  resulting  roots 
of  (41)  art  ♦  j  810  rad/sec  and  ♦  j  749  rad/sec,  where  the  J  denotes 
an  Imaginary  matter.  These  values  are  within  18  to  20  percent,  respectively, 
of  the  rocket  spin  rate  of  987  rad/sec. 

That  the  whirling  rate  during  the  spin-up  phase  Is  less  than  the  rocket 
spin  rate  has  been  documented  previously,  e.g.,  see  Figure  22  of  [1]. 

However,  no  explanation  for  the  discrepancy  was  given.  Examination  of 
this  data  shows  that  the  whirling  rate  Is  about  406  rad/sec  while  the 
spin  rate  Is  about  510  rad/sec.  This  ratio  of  whirling  rate  to  spin 
rate  of  .796  agrees  very  well  with  theoretical  results  obtained  In  this 
study.  In  fact,  the  average  value  of  the  two  whirling  rates,  810  rad/ 
sec  and  749  rad/sec,  ratloed  to  the  spin  rate  of  987  rad/sec  gives  .794. 

Other  range  data  Indicates  that,  when  the  spin  motors  are  exhausted, 
the  rocket  settles  Into  a  state  of  synchronous  whirl.  Strain  gage  data 
from  hard  mounted  launchers  support  this  contention. 


0.  PRESCRIBED  MOTION  ANALYSIS  OF  THE  PAYLOAD  SECTION 


This  section  of  the  report  will  be  an  analysis  of  a  more  prescribed 

"In  tube”  rocket  motion.  It  will  also  take  Into  account  the  actual  front 

and  rear  bearing  clearances  of  ANSSR  II.  The  basic  bent.  In-tube  rocket 

shape  Is  shown  In  Figure  4  where  f-j  and  «2  are  the  rear  and  front  bearing 

clearances  respectively.  As  shown  In  Figure  4,  a  Is  the  angle  between  the 

unbent  rocket  axis,  a  Is  the  angle  between  the  longitudinal  axis  of  the 

launch  tube  and  the  bent  rocket  axis.  Nominal  values  of  and  *2  can 

found  to  be  0.003155  Inches  and  0.002075  Inches  respectively  [5].  If  L, 

the  length  between  bearing  centers,  Is  2.1625  ft  [5]  then  the  angle  a  can 

-4 

be  found  and  Is  2.0154  x  10  rad. 

Experimental  strain  gauge  and  optical  lever  data  [1  or  2]  strongly 
Indicates  that  the  In-tube  motion  Is  one  of  synchronous  whirl,  l.e. ,  the 
rocket  spins  about  its  bent  axis  with  a  rate  of  987  rad/sec  and  It  whirls 
around  the  launch  tubes  longitudinal  axis  at  the  rate  of  987  rad/sec 
simultaneously.  This  means  that  there  Is  no  flexing  of  the  motor  section 
and  thus  no  shear  forces  due  to  flexure.  Thus,  a  somewhat  less  complex 
model  can  be  used  than  the  Timoshenko  beam  [8]  model  that  Is  needed  to 
account  for  shear  due  to  flexure. 

So,  oven  though  the  general  motion  analysis  of  the  payload  section 
yielded  a  whirling  rate  within  18  to  20  percent  of  the  actual  rocket 
spin  rate,  this  portion  of  the  report,  due  to  the  strong  experimental 
evidence  previously  cited,  will  assume  the  rocket  to  have  a  whirling 
frequency  equal  to  the  spin  frequency,  l.e.,  a  synchronous  whirling  situa¬ 
tion  prevails. 
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Bent  Rocket  Axis 
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Figure  4:  Basic.  Bent,  In-Tube  Rocket  Shape 


One*  again,  the  payload  section  Is  assuned  to  be  rigid.  This 
analysis,  due  to  symmetry.  Is  an  Instantaneous  one  and  the  particular 
configuration  In  space  Is  depicted  by  Figure  5. 

The  ground  and  body  reference  definitions  remain  the  sane  as  In  the 
previous  section  (noting  the  addition  of  the  front  bearing  clearance. 
Figure  5). 

The  Euler  angles  definitions  also  remain  the  same  but  are  placed 
under  the  constraints  that  #  *  *  9  ■  0  and  5  Is  a  constant. 

The  relation  between  the  coordinates  of  the  body  frame  and  those  of 
the  ground  reference  can  be  found  by  using  equations  (12),  (13),  (14) 
and  (15)  subject  to  the  above  constraints. 

The  terms  w,  wR,  and  a  are  as  previously  defined  and  their  respective 
values  are 


C9  (^2  +  $d 

(42) 

(n  +  C0u)e2  +  Sdaie3 

(43) 

ne2 

(44) 

The  angular  momentum  vector  of  the  payload  section  with  respect  to 
its  mass  center  is  of  the  form: 

Fcmp  3  !yy  {n  +  C9u>)*2  +  lzz  S$lj*3  (45) 

where 


w  »  the  magnitude  of  the  whirling  rate  vector,  987  > ad/sec. 
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The  sum  of  all  external  moments  about  the  payload  sections  mass  center 
can  be  found  by  substituting  equations  (42)  and  (45)  Into  equation  (20). 

The  resulting  expression  for  the  moment  vector  Is 

*oop  *  S*C#(I22  -  Iyy)  -  Ijy  -  8  S#]  (46) 

As  In  the  general  motion  analysis,  a  constant  spin  rate  of  987  rad/ 
sec  will  be  used,  second  order  terms  are  Ignored,  and  small  angles  are 
assumed. 

The  external  forces  that  create  a  moment  about  the  mass  center  of  the 

payload  section  are  shown  In  Figure  5.  These  forces  can  be  found  with  the 
use  of  equation  (23).  The  definitions  of  mp,  7,  and  r  will  remain  the 
same  but  the  expression  for  a"  and  r  will  be  written  as 

J*»x  (w  x  r)  (47) 

and 

r  *  (*  -  S#«2)  eg  +  C«f2  ®3  (48) 

Note  that 

i  *  the  linear  distance  between  the  point  "anp"  and  point  "P" 

(*  ft). 

The  resulting  expression  for  a  can  be  found  by  substituting  equations  (42) 
and  (48)  Into  equation  (47), 

a  =  w2  (S 6*2  -  S2*  *)e2  -  u2(C0«2  -  S6C#i)e3  (49) 

It  has  been  demonstrated  that  the  thrust  of  the  rocket  engine  has 
little  effect  on  the  problem  so  this  section  of  the  report  will  not 
include  it. 


The  external  moments  acting  on  the  payload  section  are  also  shown  In 
Figure  5  and  can  be  evaluated  by  substituting  equations  (23),  (46),  (48) 
and  (49)  Into  equation  (28).  Thus,  the  following  expression  Is  found  for 
the  reaction  moments  at  point  "P", 

■  -  pi  t  «2  <2  +  “2  (lyy  "  *22  ”  "p  *2)9l*j  (50) 

An  equal,  but  opposite,  moment  must  act  on  the  motor  section  and  since 
at  this  Instant  the  x  axis  and  X  axis  are  parallel,  this  moment  Is  of 
the  following  scalar  form: 


Mx  -  M  i  *2  «2  ♦  «2  (I^  -  lu‘”p  t2) 

By  aMgnlng  7  and  X,  and  by  noting  that  a  (Figure  4)  Is  equal  to 
then  equations  (9a)  and  (51)  can  be  combined  to  form, 

KO  ■  H  t  M2  +  ^(iyy  "  l2Z  ■  m^  l2)^ 

If  a  dynamic  stiffness  is  defined  as 

k  -  “2  (I>»  ■ '« -  % ,2) 
and  It  Is  observed  that 


(51) 

a 

73 

(52) 


(53) 


#  ■  a  -  a 

then  the  foilowlnq  expression  can  be  written: 
m  i  «2  c  |c* 

a  ■  - - -  i 

K  -  k 


(54) 


(55) 


This  equation  results  In  a  direct  computation  for  n.  If  a  rocket 
length  of  48.24  Inches  Is  used  and  the  rest  of  the  needed  physical  char 
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acfierlstlcs  art  obtained  free  Tables  I  and  II  then  the  resulting  value  of  a 
is  1.123  aillirad.  A  value  of  9  can  now  be  found  fro*  equation  (54)  and  Is 
0.922  ■llllrad.  Once  It  Is  realized  that  #  Is  the  "In- tube"  pitch  angle  of 

the  payload  section.  It  is  seen  that  the  presented  value  Is  In  excellent 
agreement  with  the  expert  awn tel  data  presented  In  [1  or  2]. 
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E.  HASS  CENTER  AND  PLAI  OF  THE  BENT  ROCXET 

This  section  of  the  report  deals  with  the  calculation  of  the  Mass 
center  and  the  PLAI  (Principal  Linear  Axis  of  Inertia)  of  the  bent  rocket. 
Both  will  be  found  via  a  limped  parameter  model. 

The  analysis  will  be  divided  Into  two  parts:  first,  the  mass  center 
of  the  bent  rocket  will  be  found,  and  then  the  PLAI  will  be  found  with 
respect  to  the  launcher  axis.  In  both  parts,  the  physical  characteristics 
of  the  rocket  will  be  lumped  In  four  sections:  the  motor  assembly,  the 
first  and  second  halves  of  the  shell  section,  and  the  payload  section  as 
In  Figure  6.  The  resulting  PLAI  will  give  a  mallaunch  rate  that  Is  about 
halfway  between  the  mallaunch  rates  of  200  ml  111  rad/sec  and  300  mini  rad/ 
sec  reported  In  [1]  and  [2]. 

The  mass  of  the  entire  shell  section  can  be  found  using  a  shell 
length  (L$)  of  1.871  ft.,  a  shell  thickness  of  0.038  In.,  an  Inside  shell 
radius  of  2.0  In.,  mild  steel  shell  (density  equal  to  15.23  slug/ft3), 
a  fuel  density  of  3.28  sluq/ft3,  an  inside  fuel  radius  of  1.0  In.,  and  a 
fuel  thickness  of  1.0  in  [9],  The  resulting  shell  and  fuel  mass  (m$)  Is 
0.4992  slug.  As  the  entire  mass  of  the  motor  section  is  0.7689  slug. 

Table  II, the  resulting  motor  assembly  mass  (mm)  Is  0.2697  slug.  Note  that 
the  motor  section  Is  made  up  of  the  motor  assembly  and  both  shell  sections. 

The  distance  of  the  mass  center  of  the  motor  assembly  (oi^)  from  the 
rear  of  the  rocket  (d),(see  Figure  7),  can  be  found  using  the  following 
equation, 

m  L_  +  m  (L  -  Lm  -  1/4  Le) 
d  «  m  on  s  cm  m _ s_ 

m_  v  ' 
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where 

L  •  the  distance  of  the  mss  center  of  the  no  tor  section  from  the 

CM 

rear  of  the  rocket  (0.9133  ft). 

The  resulting  numerical  value  of  d  Is  0.2556  ft. 

Figure  8  shows  the  geometry  of  the  four  rocket  sections  In  their  bent 
configuration.  The  following  symbols  will  be  defined  to  clarify  Figure  8, 

Zj  *  the  value  of  the  1  ■  1  deflection  equation  evaluated  at 

y  ■  0.333  ft.  The  numerical  value  of  Zj  Is  2.5313  x  10“4ft. 

Zj  ■  the  value  of  the  1  »  2  deflection  equation  evaluated  at 

y  *  1.269  ft.  The  numerical  value  of  Zg  Is  6.5677  x  10”4ft. 

Sa  ■  the  angle  between  the  longitudinal  axis  of  the  motor  assembly 
and  the  unbent  axis  of  the  rocket  (0.795  mllllrad). 

»  the  angle  between  the  front  half  of  the  longitudinal  axis  of  the 
shell  section  and  the  unbent  axis  of  the  rocket  (0.4312  mllllrad). 
®s2  *  *9^®  tetween  the  second  half  of  the  longitudinal  axis  of  the 

shell  section  and  the  unbent  axis  of  the  rocket  (0.7021  nrilllrad). 
Bp  *  a  *  the  angle  between  the  longitudinal  axis  of  the  pqyload  section 

and  the  rocket's  unbent  axis  (1.1236  mllllrad). 

Figure  9  shows  the  locations  of  the  mass  centers  of  the  four  rocket 
section  with  respect  to  the  (X,Y,Z)  coordinate  system.  The  mass  center  of 
the  motor  assembly  Is  described  by  coordinate  pairs  (y^,*^)*  Each  remain¬ 
ing  section  has  Its  respective  coordinate  pairs  of  (y$i»zsj)»  (yS2*zs2^* 

and  (y_,z,J.  The  numerical  values  of  each  of  these  pairs  are: 

P  P 

ym  -  0.2557  ft.  2^  -  1.942  x  10“4ft 

y$1  =*  0.8010  ft.  z$1  ■  4.551  x  10_4ft. 
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y$2  -  1.7365  ft.  zj2  ■  3.2838  x  10"4ft. 

y„  ■  2.9834  ft.  z„  -  8.7548  x  10’4ft. 

P  P 

The  mass  center  of  the  rocket  can  now  be  found  with  the  use  of  the 
following  equations. 


(57) 


where 


(58) 


y__  ■  the  y  coordinate  of  the  mass  center  of  the  bent  rocket  In 
on 

the  (X.Y.Z)  system  (2  ft) 

-  the  z  coordinate  of  the  mass  center  of  the  bent  rocket  in  the 
cm 

(X.Y.Z)  system  (2  ft) 

M  *  m,  si,  s2,  and  p  and  represents  the  motor  assembly,  the  first 

half  of  the  shell  section,  the  second  half  of  the  shell  section, 

and  the  payload  section  respectively. 

y  *  the  y  coordinate  of  the  nth  section's  mass  center  in  the  (X,Y,Z) 
n 

system  (,  ft) 

z_  =  the  z  coordinate  of  the  nth  section's  mass  center  in  the  (X.Y.Z) 
n 

system  (!?  ft) 

mn  =  the  mass  of  the  nth  section  slug). 

The  resulting  numerical  values  for  yCT  and  z^  upon  evaluation  of 
equations  (57)  and  (58)  on  1.951  ft.  and  -2.781  x  10"4  ft.  respectively. 
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The  coordinate  system  (X^.Xg.X^)  Is  shown  In  Figure  9.  The  relation-  n 

ship  between  the  coordinates  of  the  (X^.X^.X^)  system  and  those  of  the  (X»Y,Z) 
system  Is,  [j 


(59) 


r. 

!.: 

i  i 


where  the  matrix  [D]  Is  defined  as. 


>  / 

i 


CD]. 


*1  0  O' 

0  c$  -s+ 
_0  s<p  c$. 


(60) 


Note  that  $  Is  defined  by  Figure  4  and  has  a  value  of  0.2015  mllllradlans 
for  the  bearing  clearances  given  previously. 

If  X2ot  and  Xj0"  represent  the  coordinates  of  the  rocket’s  mass  center 

In  the  (X1,X2,X3)  system,  their  values,  via  equations  (59)  and  (60),  are 

1.951  ft.  and  1.1512  x  10”*  ft.  respectively. 

N  N 

The  coordinate  pairs  (X2,  X3)  can  also  be  found  with  the  use  of  equa¬ 
tions  (59)  and  (60)  and  their  values  are. 


X2  *  0.2557  ft 

X*1  •  0.8010  ft. 

Xs2  =  1.7365  ft. 
2 

XP  *  2.9834  ft. 


X3  =  2.4573  x  10”4  ft. 
X*1  -  6.1654  x  10”4  ft. 
X*2  =  6.7835  x  10"4  ft. 

X§  =  -2.742  x  10”4  ft. 


f 


f 


59 


Now  that  the  geometry  and  mass  center  locations  are  known  for  the  four 
rocket  sections  and  for  the  entire  rocket,  the  task  of  finding  the  Inertia 
tensor  of  the  rocket  can  be  undertaken.  The  Inertia  tensor  of  each  rocket 
section  will  be  found  with  respect  to  the  principal  axis  system  that  Is  loca¬ 
ted  at  Its  own  mass  center,  e.g.,  see  Figure  10.  Then  the  Inertia  tensor  of 
each  section  Is  found  In  a  coordinate  system  at  Its  mass  center,  that  Is  par¬ 
allel  to  the  (X^.X^.X^)  system.  Next,  the  entire  Inertia  tensor  of  the 
rocket  will  be  found  at  a  coordinate  system  located  at  the  rocket's  mass 
center  and  parallel  to  (X^.X^.X^).  The  final  step  Is  to  diagonalize 
this  Inertia  tensor  and  thus  find  the  bent  rocket's  PLAI. 

The  longitudinal  mass  moment  of  Inertia  of  the  first  half  of  the  shell 
section  with  respect  to  Its  principal  axis  system  can  be  found  from  the 
following  equation, 

s-jiy'y'  ■  1/2(M| r20  -  m2  r^)  +  1/2  (mn  r*fQ  -  mf2  r2fl)  (61) 
where, 

rQ  *  the  outside  shell  radius  (=  ft), 

r-j  *  the  inside  shell  radius  (f?  ft), 

r^Q  *  the  outside  fuel  radius  ft), 

rfl  *  the  Inside  fuel  radius  (j!  ft), 

and  the  quantities  m^ ,  m2,  nyj,  and  m^  are  given  by  the  following  equations, 
m-j  3  1/2  p  j  L$  r2Q  (62) 

ntj  *  1/2  p  v  Ls  r2i  (63) 

"fl  ■  1/2  »f*  Ls  r2fo 


(64) 
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"f2  “  1/2  °f  *  Ls  r2f1  (®5) 

where 

p  ■  the  mass  density  of  mild  steel  (0  slug/ft3), 

■  the  mass  density  of  the  fuel  (2  slug/ft3). 

The  transverse  mass  monent  of  Inertia  of  the  first  half  of  the  shell  section 
with  respect  to  Its  principal  axes  can  be  found  from  the  following  equation, 

s1lxV  *  sl^'z'  *  1/12  C3(,V2o  '  m2rV  +  (1/2  Ls)2(ml  ’  m2)] 

+  1/12  [3(mflr2f0  -  mf2r2fQ)  +  (1/2  Ls)2<mfl  -  mf2)]  (66) 

The  numerical  values  of  si Iy *y •  and  sl^x'x*  are  *-8731  x  10  slug/ft2 
and  0.0206  slug/ft2  respectively.  As  the  physical  dimensions  of  the  first 
and  second  halves  of  the  shell  section  are  the  same,  their  mass  moments  of 
Inertia  are  the  same. 

The  principal  axis  system  of  the  motor  section  (x',y’,z’)  is  shown 
In  Figure  11.  The  values  of  d^,  and  d$£  are  0.6577  ft.,  0.1123  ft., 
and  0.8231  ft.  respectively.  The  values  of  the  principal  longitudinal, 
and  transverse  (I^.y,  and  Ix'x'  a  Iziz>)  mass  moments  of  inertia  of  the 
motor  section  with  respect  to  (x',y',z')  are  0.01837  slug/ft2  and  0.43900 
slug  ft2  (Table II)  respectively. 

The  longitudinal  mass  moment  of  inertia  (mIy»y' )  motor  assembly 

with  respect  to  its  own  principal  axis  system,  m(x',y',z'),  can  be  found 
from  the  following  equation, 

mVy*  3  *y'y’  "  si  *y'y'  ’s2*y’y' 


(67) 
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Also,  the  transverse  «ass  moment  of  Inertle  (BIX«X«  ■  »otor 

essembly  with  respect  to  the  same  $y seem  Is  expressed  by, 

iJx'x'  "  m^z'z*  "  ^x'x*  "  sl^x'x*  ’  ^sl  ^sl  ^2  “S2^x'x'  ”  *s2^s2^2 


•  W2  '  (68) 

The  resulting  values  of  mIy'y«  mIxix»  are  8.6264  x  10”  3  slug/ft2  and 
0.1087  slug/ft2  respectively. 

The  Inertia  characteristics  of  the  payload  section  are  In  Table  2. 

The  values  of  Iy Iy ,  and  pIx.x»  are  0.0108  slug/ft1  and  0.20423  slug/ft2 
respectively. 

The  coordinate  systems  N(x.|,x2,x3)  and  an(xi*x2»x3^  are  shown  In 
Figure  12.  The  .j(xj,x2,x3)  system  Is  at  the  mass  center  of  the  Nth 
section  and  parallel  to  the  (x-j.Xg.Xj)  system.  The  ^(jc^ ,x2 *x3)  system  Is 
located  at  the  mass  center  of  the  bent  rocket  and  Is  also  parallel  to  the 
(x1,x2,x3)  system. 


Figures  13,  14,  15  and  16  show  the  geometric  relationship  between 
the  N(x',y',z’)  systems  and  the  N(x^,x2,x3)  systems. 

The  Inertia  tensor  of  each  section  can  be  found  In  Its  respective 
{x^.xj.x^  system  by  the  transformation  properties  of  the  Inertia  terms 
(See  Section  16.5  of  [10]). 

The  resulting  inertia  tensor  for  each  section  with  respect  to  Its  own 
(xj,x2>x3)  system  Is  for  the  following  motor  assembly; 


1*1  2 
m  x.jX.j  m  XjX3  *  0.1087  slug  ft 

I  -  9.6265  x  10”3  slug  ft2 

m  x2x2 

I  -  I 

m  x,x2  m  x,jX3  *  0 


1c  Re 
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_I„  ,  ■  slug  fV 

m  «• 

I  ■  8.6265  x  10"3  slug  ft3 

"  XaXa 


■  X1X2  m  X1X3 

I  -  9.6223  x  l(f5  slug  ft2 
w  ^2^3 


for  the  first  half  of  the  shell  section 


si1,,*,  ■  sl'x^  ■  °-02064  S,U9  ft 


sl'xjXj  •  4.8731  x  10"3  slug  ft2 


s,'xlx2  *  sl!x,x3  ■  0 


slIx2x3  »  9.98  x  10-6  slug  ft2 


for  the  second  half  of  the  shell  section 


s2Ix1x1  *  s2Ix3x3  3  °-02064  sl'J9  ft 


s2Ix1x2  3  s2Ix]x3  -  0 


s2!x2x3  -  -3.9575  x  10"9  slug  ft2 


for  the  payload  section 


.1  y  •  „Iy  y  -  0.20423  slug  ft" 
p  x1x1  p  x3x3 


pI  »  0.010833  slug  ftfc 

I  -I  •  0 

p  XjX2  p1x1x3 

„  -  -1.783  x  10"4  slug  ft2 
p  x2x3 

The  total  Inertia  tensor  of  the  bent  rocket  with  respect  to  Its 
an(x^,x2,x3)  system  can  be  found  with  the  following  equations, 

„  ♦  elIy  „  +  .,L  „  +  Jx  x  ♦  B »m(x?n  -  Jt) 
an  *2*2  ^  ^2^2  ^  ^2^2  ^'2*'2  P  ^2*2  *  4  » 

(69) 

+  MO?  “  x3  ^  +  ms2^x3  “  x3  ^  +  mp^x3  “  *3) 

I  *  I  .  ,  I  +  -I  +  I 
cm  x3x3  n  x3x3  si  x3x3  s2  x3x3  p  x3x3 

+  '  x2^  +  msl^x2  "  x2  ^ 

+  ms2(x°"-  x*)2  ♦  mp(x^-  xP)2  (70) 

crn^x^  "  m^x-i  +  sl^x^  +  s2*x-|Xj  +  p*x^ 

+  mm[(x™  -  x*)2  ♦  (x®-  x*)2] 

.  +  mp[(x“  -  xP)2  +  (x™  -  xP)2] 

♦  nL,[(x«  -  x3)2  ♦  (X«  -  x*2 )2 


I 

I 
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As  (xj"  -  x")*«  (x™  -  x!J)2,  then  mlx^ 


an  x^x^  an  XjX3 


(72) 


an^XgXj  "  m^gXj  +  sl^Xj  +  s2*x2x3  +  mm  ^  *3  +  "si  r2  r3 
+  ms2  r22  r32  f  mp  *2  i 


(73) 


N  M 

where  the  quantities  and  r'3  are  defined  In  the  following  equations 


rN  ■  r«  J  ♦  rf‘  k 


(74) 


and  rN  Is  the  position  vector  of  the  mass  center  of  the  Nth  section  with 
respect  to  CT(x^ ,x2,x3) ,  [10],  The  resulting  numerical  values  are: 

Iv  v  -  0.02921  slug  ft2 
cm  XjXj 


an^XjXj  3  anJx3x3 


2.2949  slug  fV 


I  w  3  -6.1663  x  10"^  slug  ft2 
an  X2*3  3 


The  angle  between  the  PLAI  axis  and  the  Xg  axis.  Figure  17,  is  given 


by 


-  2 _ I 


cm  XjX-j 


™(2W=  I  -----7 


(75) 


cm  XgX3  cm  x,x3 


The  resulting  value  for  Bp^j  is  -0.2722  millirad. 

The  instantaneous  mallaunch  rate  that  corresponds  to  this  particular 


geometry  can  be  found  from 


Geometric  Relationship  Between  the  PLAI  Axis 
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^  -  o  sin  (epLAI)  (76) 

The  resulting  value  of  Is  -  267  mllllrad/sec.  As  mentioned  previously 
this  Is  about  halfray  between  the  values  of  200  mllllrad/sec  and  300  mllll¬ 
rad/sec  presented  In  [1]  and  [2]  that  are  predicted  to  cause  the  down-range 
motion  of  ANSSR  II. 


0.3333  ft 


1.829  ft 


2. 2042  ft 


EI1 

EI2 

•  EI3 

1,006,944  1b  ft2 

208,333  lb  ft2 

1,041,667  1b  ft2 

TA8LE  I.  Physical  Characteristics  of  the  Motor  Section 
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Station  0.0  at  Nose  Tip  TABLE  II.  Inertial  Characteristics  of  the  ANSSR  II  Rocket 
Inertia  about  Center-of-Mass 
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SECTION  III 


MAL LAUNCH  AND  MALAIH  DUE  TO 
LAUNCHER  MOTION 

ABSTRACT 

Dynamical,  equations  o{  motion  one  c ieveloped  {ok  a 
multitube.  launches.  assembly  {ok  a  t pin-itabilized  Kocket. 
A  panametAic  study  is  pen{onmed  to  minimize,  the  mallaunch 
Kates  a{{ecting  missile  {light.  The  Kesults  thorn  that 
pKudent  &  election  o{  tpeci{ic  stnuctUKal  paKametens  can 
minimize  mallaunch.  Katet. 


LIST  OF  SYMBOLS 


A 

a 

B 

c 

d 


Constant  for  the  bearing  amplitude 

Average  acceleration  of  the  rocket 

Constant  for  the  bearing  amplitude 

Lateral  Stiffness  of  the  supports 

Translation  displacement  of  rocket  along  the  y-axls 
after  boost  Ignition 

Vertical  stiffness  forces 
Lateral  stiffness  forces 
Forces  due  to  pitch 
Forces  due  to  yaw 

Reversed  Inertia  force  due  to  the  vertical  acceleration 
of  launch  assembly 

Reversed  Inertia  force  due  to  the  lateral  acceleration 
of  the  launch  assembly 

Pitch  moment  of  Inertia  of  launch  assembly 

Roll  moment  of  Inertia  of  launch  assembly 

Yaw  moment  of  inertia  of  launch  assembly 

Vertical  stiffness  of  supports 

Distance  from  aft  end  of  launcher  to  rear  support 

Distance  from  rear  support  to  front  support 

Distance  from  aft  end  of  launcher  to  the  center  of  mass 
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79 


■ 

P„ 


*F 

t 

x 

xc 

y 

z 

9 

* 

+ 

u 


Center  to  center  distance  between  rocket  bearings 

Vertical  distance  from  top  of  launcher  to  the  center 
of  mass 

Mass  of  the  launch  assembly 
Magnitude  of  the  rear  bearing  force 

Magnitude  of  front  bearing  force 

Rear  bearing  force  vector 

Front  bearing  force  vector 

Time  after  boost 

Displacement  In  the  x-dl recti  on 

Distance  to  the  bore  axis  of  the  launch  tube  from  the 
C.M  In  the  x-dlrectlon 

Displacement  In  the  y-dlrectlon  of  launch  assembly 

Displacement  In  the  z-dlrectlon  of  launch  assembly 

Pitch 

Roll 

Yaw 

Operational  spin  rate  of  the  rocket 


A.  INTRODUCTION 


The  purpose  of  this  section  Is  to  present  the  results  of  an  analytical 
studty  of  the  launch  dynamics  of  a  multitube  launching  assembly  for  an 
eerodynamlcally  neutral  spin-stabilized  rocket. 

The  system  considered  In  this  section  has  two  component  parts;  the 
rocket  and  the  multitube  launching  assembly.  The  launching  assembly  Is 
composed  of  Individual  tubes  that  are  rigidly  mounted  together  to  form  a 
cluster.  Each  tube  contains  four  guide  rails  one  each  on  the  sides  of 
the  tube  and  one  each  on  the  top  and  bottom  of  the  tube  as  seen  In  Figure 
1.  The  stability  of  the  rocket  is  obtained  through  gyroscopic  effects  due 
to  the  spinning  of  the  rocket  within  bearings  mounted  integrally  with 
the  rocket  case.  This  spinning  motion  Is  attained  before  boost  (i.e.. 
Ignition  of  the  thrust  motor)  by  four  spin  motors  that  are  mounted  on 
the  aft  end  of  the  rocket.  After  the  rocket  is  brought  up  to  the  opera¬ 
tional  spin  rate,  the  thrust  motor  ignites  propelling  the  rocket  out  the 
tube. 

If  the  rocket  is  rigid  (i.e.,  the  rocket  case  does  not  bend)  and  the 
rocket's  center  of  mass  is  to  be  located  on  the  spin  axis,  then  the  rocxet 
is  theoretically  balanced.  If  either  of  these  conditions  are  not  attained, 
the  rocket  will  be  in  a  state  of  unbalance.  For  the  system  considered, 
the  unbalance  effects  are  transmitted  to  the  tube  guide  rails  via  the 
rocket's  bearings.  When  the  thrust  motor  is  ignited,  the  rocket  will 
traverse  the  length  of  the  tube  thus  causing  the  unbalance  forces  to  trans¬ 
late  along  the  axis  of  the  launch  tube.  Due  to  the  spinning  motion  of  the 


tactical  Launch  Tube 
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rocket,  the  unbalance  forces  rotate  In  a  circular  path  within  the  launch 
tube.  Thus,  the  forces  applied  to  the  launching  system  by  the  rocket 
are  equlvilent  to  a  set  of  spiraling  forces  out  the  launch  tube.  This 
type  of  Motion  by  the  rocket  within  the  tube  can  cause  the  launching 
assembly  to  oscillate.  If  the  oscillations  are  severe  enough,  they  can 
Impart  to  the  rocket  a  malalm  and  a  mal launch  rate  at  the  end  of  guid¬ 
ance. 

The  goal  of  this  Investigation  is  to  determine  the  natural  frequencies 

of  the  launching  system  and  to  determine  If  the  assembly  has  any  undue 

rotational  motions  at  end  of  guidance  severe  enough  to  disturb  missile 

launch.  The  theoretical  work  has  been  directed  toward  formulation  and 

solution  of  the  dynamical  equations  of  motion  for  the  multitube  launching 

assembly.  The  rocket  and  the  assembly  are  treated  as  discrete  parts  with 

the  rocket  contributing  only  to  the  motion  of  the  assembly  by  the  amplitude 

and  the  spin  rate  of  its  bearing  reactions.  In  an  earlier  study  It  was 

found  that  the  points  of  application  of  the  force  vectors  of  the  front 
and  rear  bearings  are  180  degrees  out  of  phase  (see  Section  II).  In  the 

mathematical  expressions,  the  system  considered  Is  one  with  two  supports 

equally  spaced  longitudinally  from  the  center  of  mass  as  seen  in  Figure  2. 

Each  support  has  a  given  lateral  and  vertical  stiffness.  As  the  two 

supports  are  equally  spaced  from  the  center  of  mass  of  the  launcher 

assembly,  the  dynamical  equations  of  motion  are  decoupled.  As  a  result 

of  this  decoupling,  mathematical  expressions  can  be  developed  in  a  straight 

forward  manner  to  determine  the  theoretical  position  and  state  of  motion 

of  the  cluster  assembly  during  launch. 

From  the  mathematical  expressions  for  the  natural  frequency,  the 
position,  and  the  motion  of  the  assembly,  parametric  results  are  obtained 


Cluster  Assembly  with  Supports 
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and  are  presented  in  part  D  of  this  section.  Thus,  given  certain  para¬ 
meters  that  affect  the  stability  of  the  assembly,  notably  the  acceleration 
of  the  rocket,  the  amplitude  of  the  bearing  forces  the  spin  rate  of  the 
rocket,  and  the  stiffness  of  the  mounting  hardware,  the  malaim  and  the 
mal launch  rates  at  the  end  of  guidance  can  be  determined  and  accounted 


for. 


B.  MATHEMATICAL  FORMULATION  OF  DYNAMICAL  EQUATIONS  OF  MOTION 


In  developing  the  differential  equations*  which  describe  the 
Motion  of  the  cluster  assembly,  several  physical  dimensions  are  used. 
For  a  complete  list  as  to  all  of  these  physical  dimensions,  reference 
Is  made  to  Table  I.  It  Is  assumed  that  the  Individual  tubes  of  the 
cluster  assembly  are  rigid  with  respect  to  one  another.  The  center 
of  the  mass  Is  considered  fixed  and  the  supports  are  equally  spaced 
with  respect  to  the  center  of  mass. 

The  coordinate  system  to  which  all  motion  and  position  Is 
referred  Is  a  right-handed  cartesian  coordinate  system,  X,Y,Z,  with 
the  origin  0  of  the  system  fixed  at  the  center  of  mass  In  such  a 
position  that  the  Y-axis  lies  parallel  to  the  bore  lines  of  the  launch 
tubes . 

From  the  free  body  diagrams.  Figures  3  and  4,  the  summation  of 
forces  In  the  negative  z-dl recti  on  gives 

£F_Z  «  0  »  -  %  +  Fi  ♦  Fi  -  Rp  ♦  F2  (1) 

and  the  summation  of  forces  In  the  positive  x-dlrectlon  gives 

£F„  ■  0  -  -  RB  *  F3  ♦  Fj  -  Kf  *  F,  (2) 

where  Fj  and  Fj  are  the  reversed  Inertia  forces  due  to  the  vertical  and 
lateral  stiffness  of  the  launcher  assembly  such  that 
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C.M.  I 

3 - «fFJ 

U — H 

Free  Body  Diagram  Number  2 
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TABLE  I 

Physical  Dimensions  of  Launch  Tube 

12  Support  Attachment  Spacing  14.0  In 

13  Transverse  Distance  to  Center  24.495  In 

J  of  Mass  Forward  of  AFT  End  of 

Launch  Tube 

14  Rocket  Bearing  Spacing  25.690  In 

I  Pitch  Moment  of  Inertia  13.6  slug  ft2 

I  Yaw  Moment  of  Inertia  13.6  slug  ft2 

m  Mass  of  Launch  Assembly  5.49  slugs 

a  Average  Acceleration  of  3122.316  ft/sec2 

Rocket 

t  Time  after  Boost  Ignition 

u  Rocket  Operational  Spin  Rate  158  cps 


88 


+  o3  -  d)  irB  -  o3  -  d  - 14)  irF 

From  the  geometry  of  the  system  and  the  stiffnesses  of  the 
supports,  the  following  physical  relationships  are  obtained: 


F-|  *  Fg  *  kz 

(6a) 

F3  =  F4  »  F5  -  F6  -  cx 

(6b) 

90  . 

F7  ■  k  (z  -  h  12  sin  $  ) 

(6c) 

F8  ■  k  (z  +  Jj  12  sin  0  ) 

(6d) 

Fg  *  c  (x  -  h  12  sin 

(6e) 

F'jo  ■  c  {a  :  \  '2  s'n  ^ ) 

{$*) 

for  small  angles  of  rotation. 

F7  ■  k  (2  -  h  12  0  ) 

(6g) 

Fg  *  k  <z  +  H  12  $  ) 

(6h) 

Fg  *  c  (x  -  %  1 2  0  ) 

(6j) 

F10  •  c  (x  +  H  12  iff  ) 

(6k) 

The  bearing  reactions  are  cyclic  and  revolve  within  the  launch 

tube  at  a  rate  approximately  equal  to  the  spin  rate  of  the  rocket  [3]. 

The  front  and  rear  bearing  reactors  are  180°  out  of  phase.  Thus, 

the  bearing  reactions  can  be  represented  by  the  sinusoid  of  functions 

which  results  in 

TTp  *  -  Pp  sin  u  t  (7a) 

fig  =  Pg  sin  cj  t 

(7b) 

where  Pp  and  Pg  are  the  magnitudes  of  the  front  and  rear  bearing 
reactions  respectively  and  wis  the  spin  rate  of  the  rocket. 

The  distance,  d,  that  the  rocket  has  moved  out  the  launch  tube 


is  giver  by  the  equation 
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d  -  h  at2  (8) 

where  a  Is  the  average  acceleration  of  the  rocket  and  t  Is  the  elapsed 
tine  after  boost. 

Fron  Newton's  second  law  of  motion,  the  Inertia  forces  asso¬ 
ciated  with  the  cluster  assembly  due  to  acceleration  of  the  center  of 
nass  vertically  and  laterally  are 

Fj  ■  ui  (9a) 

Fj  ■  mx  (9b) 

Equations  (6a,  6b,  6g-k,  7a,  7b,  8,  9a,  3b)  substituted  Into 
equations  (1-5)  result  In  the  dynamical  equations  of  motion. 


mz  +  2kz  -  (PB  -  PF)  sin  wt  (10) 

mx  +  2kz  ■  (PB  -  Pp)  sin  ut  (11) 

Iyy^  +  2c1  (jX  «  (Pg  +  PF)(zc  coswt 

+  Xc  slnwt)  (12) 

I***  ♦  Jj  fyO  *  -C(l3  -  14)  PF  ♦  13  PB] 

+  slnwt  +  Jj  (Pg  -  PF)  at2  slnwt  (13) 

+  slnwt  +  *s(PB  +  PF)  at2  slnwt  (14) 


Equation  (12)  has  two  dependent  variables,  #  and  x,  and  one  Independent 


C.  SOLUTIONS  OF  DYNAMICAL  EQUATIONS  OF  MOTION 


The  solution  to  the  dynamical  equations  of  motion  derived  In  the 
preceding  section  will  be  given  In  two  parts,  the  first  being  the 
natural  frequency  of  the  cluster  assembly,  which  will  be  obtained  from 
the  unforced  part  of  the  equations  of  motion.  The  second  solution  will 
be  used  to  obtain  the  position  and  motion  of  the  assembly  as  a  function 
of  time  during  launch. 


Part  1 

The  characteristic  equations  are  obtained  by  setting  the  unforced 
part  of  the  equations  of  motion  equal  to  zero  which  yields 


(ms2  +  2k)z  =  0 

(17) 

(ms2  +  2c)x  =  0 

(18) 

(Iyys2  +  2cl|)0=  0 

(19) 

uXxs2  +  0 

(20) 

d2zs2  +  Jsl|c)^=  0 

(21) 

The  natural  frequencies,  5,  are  obtained  by  a 

characteristic  equation  which  gives 

solution  of  each 

S1  *  ±  yj- 2k/m  ’ 

(22) 
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Sg  ■  *y-2c/m 

(23) 

Sj  .  2cl|/Iw 

(24) 

s4  -  ±^-lg  k/2  1^ 

(25) 

s5  -  ifl*  c/2  IZJ 

(26) 

Part  2 

Each  of  the  dynamical  equations  of  motion  Is  a  second  order 
linear  differential  equation  with  constant  coefficients.  The  method  of 
variation  of  parameters  Is  used  to  solve  the  equations.  Three  different 
solutions  are  required  due  to  differences  In  the  forcing  functions  In 
the  equations  of  motion.  Equations  (10)  and  (11)  are  of  the  same  form 
and  equations  (13)  and  (14)  are  of  the  same  form.  It  Is  possible  to 
obtain  a  general  solution  for  each  set. 

The  general  form  of  equations  (10)  and  (11)  Is 

wy  ♦  ky  ■  A  slnmt 


where  m,  k,  and  A  are  constants.  This  Is  a  second  order  linear  differ¬ 
ential  equation  with  constant  coefflents  with  a  sinusoidal  forcing 
function,  y  Is  obtained  by  the  method  of  variation  of  parameters  as 
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equation  (13)  Into  equation  (30)  as 


FT7 

»ii7 


M(13-l«)Pf  ♦  )3  PB]  12  »  !„■(>,  *  fj) 


L  k  -  2  I„«2 


(1?  k  -  2  I  «2)2 

'  2  XX 


t  32  I„2  »2  .(P„  +  PF) 
('2  k  '  2  ‘xx"2)2 


sin 


jsr. 

’  2  I 

xx 


.  *(pb*pf)  .  .  ,  » 'xx-Ce  ♦  pp> 

♦  rr~ - : - r~  t*  sln«t 


>2  ‘  *  2  ‘xx"2 


(‘1  “  *  2  ‘xx"2)2 


t  C0Sto»t 


~2[(13  -  l4)Pf  ♦  13  P8]  ,  2  ‘,x2<Pg  ♦  Pf> 

if  k  -  2  I**-2  (13  k  -  2  Ijm"2)2 


32  I„2  «*3a(P6  ♦  PF) 
('2  k  ’  2  ‘xx"2)2 


slnwt 


(31) 


^  Is  obtained  by  a  substitution  of  the  variables  and  constants  of 
equation  (14)  Into  equation  (30)  as 


T 

#Visf 


2u[(l3-l4)Pf  *  13  PB]  ,  12  tj  Iz2a(Pg  *  PF) 

(122  c  -  2  I22"2)2 


'I  c  -  2  I22"2 


+ 


32  I2Z2^3a(PB  ♦  PF) 
(If  C  -  2  Iz2w2)3 


t 


(Equation  continued  on  next  page.) 
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a(PB  +  PF)  „  8  I72wa(PB  +  Pp) 

+  - IE - -  t2  sin wt  -  —5 - — —  t  coswt 

l|  c  -  2  I22«2  (l|  c  -  2  I22«2)2 


2t(13-14)PF  +  ^3Pb3  ,  4  !22a(PQ  +  PF>< 

'  l|c-2  1„u*  *  O2  c  -  2  lzy 


32  I„2«3a(Pa+PF) 

+  — r-“ - 8— 1 —  slnwt 

0|  c  -  2  IZ2W2)3 


Equation  (16)  Is  of  the  form 

■jy  ♦  Icy  »  A  cos«t  +  B  slnwt 
The  same  method  also  serves  to  obtain  a  solution  for  y  as 


y'V^ fcr)  ,,nV^t-(j±iW*t  i 


B 

k-fflw2 


•slnwt  +  ( — ~r  )  coswt  (33) 

\k-1tiw2  / 

0  follows  by  a  substitution  of  equation  (16)  Into  equation  (33) 


K 

\  <PB  '  PF> 

1 2  c  if 

2  c  lj  -  Iyyw2 

'2  c  11 


(Equation  continued  on  next  page) 
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.r  <w_i 

2  c  l|  -  I w«2  J  I  Iw  [2  c  l|  -  I nu- 


•  [Xc  slnwt  +  Zc  coswt]  (34) 

The  stiffness  of  the  mounting  hardware  Is  considered  such  that 
the  lateral  stiffness  and  the  vertical  stiffness  are  equal.  The  equa¬ 
tions  (31  and  32  are  Identical  except  for  the  principal  moments  of 
Inertia.  If  the  moments  of  inertia  are  approximately  equal*  then 
equation  (32)  becomes  Identical  to  equation  (31).  Equations  (28  and  29) 
are  Identical  as  the  vertical  and  lateral  stiffness  are  equal.  Thus* 


substitution  of 


k  »  c 


Into  equations  (28,  29,  31,  and  32)  yields 


X  •  Z  « 


■vsIKS)  -vf.  >» 


j  0  PTT  ^(VW  *  W  . 

V  l|  k  \\  k  -  2  l^u2  (l|  k  -  2  IjujW2)2 


32  *xx2  w3a(Pg+Pp)  ]  fifk  »<!■„♦  PF)  _ 

■f  — - - -  sin*/ - t  +  -5 - t2  slnwt 

(lg  k  -  2  Ij^W  )  f  2  lg  k  -  2  Ixx«2 


12  k  -  2  IXx"2 


(Equation  continued  on  next  page.) 
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8  Ixxi»>a(PB+  PF) 

'  t'i  k  - 2  V2 


2tfl,-WW  (  4  l„a{PB^F) 

.1  k  -  2  '„"2  0*  k  •  2  >»"2)2 


+ 


32  a(PB+  PF) 

Oi k  - 2  V»2>s 


slnwt 


(36) 


For  the  time  rate  of  change  of  pitch, 9  ,  and  displacement. 


0m  0*  ^<0) 


2«C(13-14)Pf  ♦  Vb3  , 
l|  k  -  2  (1 }  k  -  2IXX«2)2 


,  32  I^3  a(PB^F)  ‘ 
(l|  k  -  2  I^w2)3 


2  a(PB+PF) 

- - i 

l|  k  -  2  I^w2 


a<W 

•  slnwt  +  — - t2  coswt  + 

i| k  - 2 


8wlxx»(  V’f) 

(l|  k  -  2  I^w2)2 

»  m 


2C(13-14)Pf+13PR)]  ^  Ixx  a(PB+Pc) 

[Wt  slnwt  -  coswt] - —  -4-  —  J  8  ♦ - ^ - 

i  k  -  2  Ixxw2  (if  k  -  2  i^2)2 


32  I  2w3  a(PB+PF)  * 

+ - — - S-J—  wcoswt 

(1|  k  -  2  I^w2)3 

m 


(37) 


and 
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Equations  (35,  36,  37,  and  38)  give  the  displacement,  malalm, 
■al launch,  and  rate  of  displacement  as  a  function  of  the  parameters, 
support  stiffness,  average  acceleration,  bearing  reactions  and  spin 
rate. 


0.  ANALYSIS  AND  CONCLUSIONS 


From  the  solution  of  the  dynamical  equations  of  motion,  a  para¬ 
metric  study  has  been  done  to  determine  the  malalm  and  mallaunch  rates 
at  end  of  guidance.  The  results  of  this  study  are  presented  in 
Figures  8  thru  21  for  the  multi  tube  launcher  parameters  as  listed  In 
Table  I. 

In  this  study  the  following  parameters  are  varied:  the  average 
acceleration,  the  bearing  reactions,  the  support  stiffnesses,  the 
operational  spin  rate  and  the  time  from  Initial  boost  until  end  of 
guidance.  The  average  acceleration  was  varied  from  1000  ft/sec2  to 
4000  ft/sec2  as  shown  In  Figures  8  and  9.  Increasing  the  average 
acceleration  (l.e..  Increasing  the  thrust)  decreases  the  time  to  end 
of  guidance,  but  Increases  the  mallaunch  rate  at  end  of  guidance. 

By  assuming  the  front  and  rear  bearing  forces  to  be  equal,  the 
displacement,  z,  and  the  rate  of  change  of  displacement,  z,  are  zero. 
Increasing  the  bearing  reactions  equally  causes  the  magnitude  of  the 
malalm  and  mallaunch  rate  to  Increase  as  shown  In  Figures  10  and  11. 

Figures  12  and  13  show  typical  solutions  of  the  dynamical 
equations  of  motion  for  support  stiffnesses,  k  3  0.0  to  k  *  ®  .  As 
the  natural  frequency  of  the  cluster  assembly  approaches  the  critical 
spin  rate,  that  Is 

2  Ixx"2 
k  — — - 
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the  magnitude  of  the  malalm  and  ptal launch  rate  approaches  Infinity 
making  the  launcher  assembly  Increasingly  unstable  In  these  areas. 

The  area  of  study  Is  confined  to 

2 

k «  - 

Figures  14  and  15  show  malalm  and  mal launch  rates  as  functions  of 
time  for  support  stiffnesses  equal  to  zero. 

The  malalm  Is  a  finite  value  for  zero  stiffness  at  the  end  of 
guidance,  but  Is  zero  before  boost.  Figure  16.  The  mallaunch  Is  also 
finite  at  end  of  guidance.  Figure  17,  but  Is  zero  at  time,  t  ■  0.0. 

After  boost,  the  mallaunch  Increases  negatively  then  becomes  positive 
at  end  of  guidance.  At  end  of  guidance  the  mallaunch  rate  tends  to 
become  Increasingly  unstable  for  large  bearing  reactions  and  zero 
support  stiffnesses. 

Equations  (36)  and  (37)  are  periodic  functions  of  the  operational 
spin  rate.  Figures  18  and  19,  hence  It  Is  possible  to  minimize 
either  malalm  or  mallaunch.  The  period  for  mallaunch  Is  approximately 
1.4  cycles.  The  operational  spin  rate  need  vary  only  0.35  cycles  per 
second  at  the  end  of  guidance  for  the  mallaunch  rate  to  go  from  zero 
to  a  maximum.  For  a  missile  with  an  operational  spin  rate  of  158  cycles 
per  second,  the  deviation  from  zero  to  maximum  mallaunch  is  0.225. 

Equations  (36)  and  (37)  are  also  periodic  functions  of  time. 

The  mallaunch  rate  and  the  malalm  become  Increasingly  large.  Figures  20 
and  21  as  the  missile  approaches  end  of  guidance.  By  close  selection 
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of  the  average  acceleration  of  the  rocket*  the  tine  at  end  of 
guidance  My  be  selected  such  that  the  nal launch  rate  Is  zero. 

As  the  mallaunch  rate  1«  the  time  rate  of  change  of  malalm*  it 
will  not  be  possible  to  null  both  variables  at  the  same  time.  The 
malalm  affects  only  the  point  of  Impact  of  a  launched  missile  and 
can  easily  be  accounted  for  with  Initial  aim.  The  mallaunch  rate 
affects  the  actual  trajectory  of  flight,  and  hence,  should  be  nulled. 

The  performance  of  the  missile  Is  very  sensitive  to  any  changes 
In  Its  operational  spin  rate.  The  bearing  forces  Increase  with  the 
square  of  the  operational  spin  rate.  The  time  to  end  of  guidance  Is 
a  direct  function  of  the  average  acceleration,  and  hence,  the  thrust 
and  the  operational  spin  rate  of  the  rocket  must  be  tightly  controlled. 


Figure  8.  Mel aim  versus  Acceleration 
{  U  ■  158  cps,  Pa  ■  50  lbf) 
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Figure  11.  Mai launch  Rate  versus  Bearing  Reaction 


Figure  12.  Maltln  versus  Stiffness  »t  End  of  Guldtnce 


Malal 


Figure  16.  Malalm  versus  Stiffness  at  End  of  Guidance 
(<j  •  158  cps,  a  «  3122.316  ft/secz) 


Mallaunch  Rate  versus  Spin  Rate 
at  End  of  Guidance  (PD  *  50  Ibf 


Mai launch  Rate  versus  Time 
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F.  APPENDIX  A  COMPUTER  PROGRAM 


C  MAIN — MALAIM  ANO  MALtAUNCH  VERSUS  STIFFNESS 

OIMENSION  PI TCHOI 1*498)*  PITCHK 1,498) *  AKlll,498). 
1PBHW49B) 

26  FORMATI4E15.4) 

ON  tS  ?H£  OPERATIONAL  SPIN  RATE  OF  THE  ROCKET.  (RAD/S=C> 
0M»158.0*2.*3. 14159 


PF  ANO  PB  ARE  THE  MAGNITUDES  OF  THE  BEARING  FURCFS 
OF  THE  ROCKET'S  FRONT  AND  REAR  BEARINGS  RESPECTl VELY. 
ILBFI 

PB*53.0 

PF-50.0 


AM  IS  THE  MASS  OF  THE  MULTITUBE  LAUNCHING  ASSEMBLY • 
(SLUGS) 

AM*223.5/32.2 

AJ  IS  THE  MOMENT  OF  INERTIA  ABOUT  THE  X-AXIS. 
(SLUGS-FT-FT) 

AJ-li.6 


L2  IS  THE  DISTANCE  BETWEEN  THE  MOUNTING  LUGS  OF  THE 
LAUNCHER  ASSEMBLY.  (FT) 

U2*14./12. 


03  IS  THE  DISTANCE  FROM  THE  AFT  END  OF  IMF  ROCKET  TO 
THE  CENTER  OF  GRAVITY.  (FT) 

03*24.493  /12. 


D4  IS  THE  CENTER  TO  CENTER  DISTANCE  BETWEEN  THE 
BEARINGS.  (FT) 

04*25. 69/12. 


C 

C 
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C  A  IS  THE  AVERAGE  ACCELERATION  OF  THE  ROCKET  AS  IT  IS 
C  PROPELLED  OUT  THE  TUBE.  (FT/SEC/SEC) 

C 

A-3122.316 

C 

00  15  K-1,1 

AK  IS  THE  STIFFNESS  OF  THE  COUNTING  LINKS.  ( LBF/FTl 
AK»1000.0 
UO  11  J* l » 498 

T  IS  THE  ELAPSED  TIME  FROM  INITIAL  BOOST,  f  SEC ) 

T«0. 7 

AB«2.*AJ/(02**2*AKI 
AC* ( 03-04 )*?F*03*PB 
AD*A*(PB*PF) 

AG»02**2*AK-2. *AJ*QM**2 

L  IS  THE  VERTICAL  DISPLACEMENT.  (FT) 

Z*-$QRT(AM/(2.*AK) )*( ( PB-PF ) *0M/ { ? . *AK-AM*0M**2 ) I  * S I N { 
l  SORT  (2.*AK/AM)*T)*(  (  PB-PF  )  /  (  2.  * AK- AM*0M**2  )  I  *S  I N  ( 0>V*T ) 

I jI  IS  THE  RATE  OF  VERTICAL  01 S°L ACFMENT .  (FT/SEC) 

02  *(PB-PF)  *0M/  (2.*AK-AM*'3M**2)*(CUS(0M*T)-CQS($jRT 
1(2. *AK/ AM ) *T ) ) 

pitch  is  the  malaim.  (rad) 

PITCH  =  SQRT(  AB)*(2.*OM*AC/AE*l?.*OM*AJ*AU/A6**?-t-32.*AJ 
l**2*CM**3* AD/A  6**3 )*S  IN(SOR  T(  1./A3  )*T)«-AO/A?*T**?*Sl  'i 
2(0M*T )-8,*AJ*aM*A0/AF**2*T*CnS(0M*T)-(2.*AC/AE*4.*AJ* 
JAO/AE**2  +  32.*A,J**2*OM**3*AO/A6**3)  *SlN(riM*T  ) 

OPITCH  IS  THE  MALLAUNCH  RATE.  (RAO/SEC) 

UP  ITCH  a(2.*OM*AC/Afc*12.*OM*Aj*AH/AE**2*32.*AJ**2* 
10M**3*A0/AE**3)  *COS(  SORT(  1 .  /AP  )  *T  )  *2 .  *AU/A  6*T*S  1 1,(  0*-*  T 
2)  ♦UM*A0*T**2/AE*C0S(0H*n*  (  8.*AD*AJ*0.^/A6**2  )*(;)•*  T* 
3SIN(0M*T)-C0S(0M*T  J  )-(  2.*AC/AF*4.*AJ<-AD/A  =  **2  +  32  .*AJ** 
42*0M**3*A0/AE**  3  )  *0***C,JS  ( UM*T  ) 

PITCH0(K,J)*0P!TCH 
PlTCHl(K,J)=piTCH 
PBKK.J  )*PB 


/ 
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AK 1 (K*J ) »AK 
AK  »AK4-2000.0 
II  CONTINUE 
P8-P84225.0 

15  CONTINUE 

00  16  K«l,l 
00  16  J»l,498 

16  WRlTEt6,26)PITCHD(K,J) .PITCHlf K,J) ,PtU(K, J)*AKl(K,J) 
CALL  SPPLOTIAKl.'STIFFNESS  ( LBF/FT > • , lfl , PI TCH1 , 

I'MALAIM  (RAO) *  1 12* 

2  *F  IGURE  8  MALA  I M  VS.  STIFFNESS  AT  EOG' » 36*  5.C,  3.  3.  C-, 
11*498) 

CALL  SPPL0TIAK1, 'STIFFNESS  I  LBF/FT »', 1 B, »I TCHD, 
I'NALLAUNCH  I  RAO/SEC )•* 19, 

2'FIGURE  14  MALLAUHCH  RATE  VS.  STIFFNESS  AT  EOG', 45, 
3S..C, 1.0,0, 1,498) 

CALL  PL0T(0. 0,0. 0,999) 

STOP 

ENU 
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SUBROUTINE  SPPLOT ( X , XHOR ,KX, Y, YHOR ,KY,FHOR  «KF , 
iALENX,ALENY,LINTVP,NROW,NCOU 


SPPLOT  IS  A  SUBROUTINE  THAT  WILL  PLOT  NROW  NUMBER  UP 
LINES  PER  GRAPH. 

THE  ARGUMENTS  OF  THE  LIST  AREO 

X  NAME  OF  THE  ARRAY  CONTAINING  THE  ABSCISSA  OR  X 
VALUES. 

Y  NAME  OF  THE  ARRAY  CONTAINING  THE  ORDINATE  OR  Y 
VALUES. 

XHUR  AND  KX  DESCRIBE  THE  DESIRED  LEGEND  FOR  Th= 
X-AXIS. 

XHOR  IS  THE  LITERAL  EXPRESSION  OF  THE  LEGEND  TO  Re 
PRINTEU. 

MX  IS  THE  INTEGER  NUMBER  OF  CHARACTERS  IN  XHOR. 
YHOR  AND  KY  DESCRIBE  THE  DESIRED  LEGEND  FOR  TH~ 
Y-AXIS. 

YHDR  IS  THE  LITERAL  EXPRESSION  OF  THE  LEGE  NC.  TO  h- 
PRINTEO. 

KY  IS  THE  INTEGER  NUMBER  DF  CHARACTERS  IN  YHOR. 
ALENX  IS  THE  DESIRED  LENGTH  OF  Th=  A-AXIS. 

ALENV  IS  THE  DESIRED  LENGTH  OF  THE  Y-AXIS. 

LIMTVP  DETERMINES  THE  TYPE  OF  LINE  YOU  GET. 

IF  LINTYP  GREATER  THAN  0,  YOU  GET  A  STRAIGHT  1 1  *|F 
CONNECTING  A  SYMBOL  PUT  DOWN  AT  EVERY  L I NTVPTM 
POINT. 

THUS.  IF  LINTYP  WERE  EQUAL  TO  2,  YOU  WOULD  GET  A 
LINE  CONNECTING  EVERY  SECOND  POINT. 

IF  LINTYP  WERE  EQUAL  TO  3.  YOU  WILL  GET  A  Lilt 
ONLY. 

IF  LINTYP  LESS  THAN  0,  YOU  WOULD  GET  A  SYYoJL 
EVERY  LINTYPTH  POINT. 


DIMENSION  X( NROW.NCOL I »  Y( MkOw.NCOL ) «  XHOR (231. 

1YHCR(20),FH0R(1C0> 

CALL  GSIZ5I11. 0,11.0, 1121) 

CALL  SYMBOL (2. 5, 1.0, . 1 5,FHCR , 0. 0,KF ) 

CALL  PLOT (2. 5, 5.3, 3) 

CALL  PLOT  ( 7.5, 5. S, 2) 

CALL  PLOT (7. 5, 2. 5, 2) 

CALL  PLOT! 2. 5, 2. 5, -3) 

CALL  SCALE  21  X ,  ALENX  ,N«Ov»,  NCOL ,  FRX  ,  DLX ) 

CALL  SCALF2(Y,ALCNY, NROw.NCOL *  FRY, PLY ) 

CALL  AX  I  S(  0.0, 0.0,  XHOR  ,-KX  ,  ALENX  » 0 . D.P3X ,DLX ) 
CALL  AXIS! G.O, D.O.YHOR, KY, ALENY,<n.O,FRY,JLY) 


i 

f 
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00  32  J* I  *  NROW 

If ( J  ,N£.  1)  CALL  PLOTIO.0,0.0,3) 
INTEQ-J 
IPEN  •  3 
I COO 6  *  -1 
NT  •! A8SIL INTYP  I 
IE  I L INTYP) 7, 6, 7 
6  NT*1 
T  NF*l 
NA  >  NT 
KK>1 

IE  ILINTYP)  11,12*13 

11  I PENA  *  3 
ICOOEA  «  -1 
LSW  *  1 

SO  TO  15 

12  NA*NCOL 

13  IPENA  *  2 
ICOOEA  *  -2 
LSW-0 

15  00  30  I *1 , NCOL 

XN»IX(J,I)-F«X»/DLX 
YN*I YI J, I )-FRY) /OLV 
IF  (NA-NT)  20,21,22 

20  IF  I LSW )  23,22,23 

21  CALL  SYMBOL  ( XN , YN,0.Q8, I N TEC ,0.0, ICOOE ) 
NA  *  1 

GO  TO  25 

22  CALL  PLOT  (XN,YN,IPEN) 

23  4A  *  NA  ♦  1 
25  NF  *  NF*KK 

ICOOE  »  ICOOEA 

30  IPcN  ■  IPfcNA 

31  CONTINUE 

32  CONTINUE 
RETURN 
ENr 
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c 

SUBROUTINE 

ARRAY 

SCALE2I ARRAY ,AXLEN, NR OW.NCJL, FIRST, DEL) 
NAME  OF  ARRAY  CONTAINING  VALUES  Jn 

c 

c 

AXLEN 

SCALED. 

LENGTH  IN  INCHES  OVER  WHICH  ARRAY  IS  TO 

c 

c 

NPTS 

BE  SCALED. 

NUMBER  OF  POINTS  TO  RE  SCALED. 

c 

INC 

INCREMENT  OF  LOCATION  OF  SUCCESSIVE 

c 

POINTS. 

DIMENSION  ARRAYINROW.NCOLI ,SAVE<7) 

SAVE  1 1 ) -l. 0 

SAVE (21*2.0 

SAVE! 31-4.7 

SAVE  141-5.  7 

SAVE  I 5 ) -B. 7 

SAVE  <61-10.0 

SAVE! 71*20. 

FAO-0.C1 

inc«i 

NPTS-NCOL 
K* I ABS ( INC) 

N«NPTS*K 
V1-1.0E60 
VF  — l.OE-67 
DC  26  JJ-l.NROW 
YC-ARRAYI JJ,1) 

VN-YO 

DO  26  l-l,N,K 
VS- ARRAY ( J  J  «  I  ) 

IF  (YO-YS)  22,22,21 

21  YO-YS 

GO  TO  25 

22  IF  (YS-YN)  25,25,24 
24  YN-YS 

il  I F ( YN  .GT.  YF ) YF-YN 
1 F  ( Y  0  .LT.  YlJYl-YG 
26  CONTINUE 
YC-Yl 
YN-YF 
FIRSTV-YO 
IF  (YO)  34,35,35 

34  FAU-FAD-l.O 

35  OELTAV»( YN-F IRSTV) /AXLCN 
IF  (UELTAV)  70,70,40 

I  * ALOGl 0  (DELTAV)>107 
P=  I0.0**(  1-1000  ) 

OELTAV-OFLTAV/P-O. 01 
UO  45  1*1,6 

IS*  I 
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IF  (SAVES ! I -DELTA VI  45*50,53 
45  CONTINUE 
50  OELTAV«SAVE(  IS)*P 

FIRSTV«06LTAV*AINT(VO/OELTAV*FACI 
T»F1RSTVMAXLEN40.0D*0ELTAV 
IF  (T-VNI  55,57,57 

55  F!R$TV»P*AINT(Y0/P+FAD) 

T-FIRSTV»(AXLEN*.Oi)*OELTAV 
IF  (T-YNI  54,57,57 

56  IS-ISM 

GO  TO  50 

57  FIRSTV-FIRSTV-AINTM AXIEN* ( F l RSTV-YN) /DFLTAV ) /2. 0 ) * 
10ELTAV 

IF  (YO*FIRSTVI  58,5ft, 59 

58  FI«STV»3.0 

59  IF  (INC)  61,61,65 

61  FIRSTV-F!»SmAINT(  AXLEN*.  5)*0£LTAV 
OELTAV—OELTAV 
65  FIRST-FIRSTV 
UEL*0ELTAV 
67  RETURN 

70  OELTAV-2.0*FIRSTV 

O£LTAV«A0S(OELTAV/AXLEN)+1. 

GO  Tu  43 
ENO 


